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ABSTRACT 


Heat-transfer  measurements  were  made  for  filmwise 
condensation  of  steam  on  smooth  tubes  and  finned  tubes  with 
rectangular-shaped  fins.  The  influence  of  fin  root  diameter  on  a 
single  horizontal  finned  tube  was  examined  by  taking  data  for  two 
families  of  tubes  whose  fin  root  diameters  were  12.7  and  19.05  mm. 
Each  family  of  tubes  had  fins  that  were  1  mm  thick  and  1  mm  high, 
but  with  fin  spacings  oi  u.25,  0.50,  1.0  1.5,  2.0,  and  4.0  mm. 


To  compare  resuits  being  obtained  in  a  concurrent  research 
project  with  Queen  Mary  College  (University  of  London),  data  was 
collected  on  two  tubes  whose  fin  root  diameters  were  12.7  mm,  and 
which  had  fin  thicknesses  oi  0.5  mm,  fin  heights  of  1.59  mm,  and  fin 
spacings  of  1.0  mm.  A  comparison  of  the  enhancement  ratios  (based 
or.  constant  vapor-side  temperature  drop)  with  past  Naval 
Postgraduate  School  (NPS)  data  showed  good  agreement. 

The  apparatus  was  then  modified  in  order  to  study  the  effects  of 
condensate  innundation.  Preliminary  flow  visualization  studies  of  the 
condensate  motion  were  made.  On  finned  tubes,  the  condensate 
falling  on  the  tubes  has  a  tendency  to  be  channeled  between  fins  and 
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does  not  spread  out  axialiy  as  occurs  on  smooth  tubes. 
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I.  INTRODUCTION 


A.  BACKGROUND 

As  technological  advancements  continue  to  enter  naval  vessels  in 
the  form  of  more  powerful  and  capable  weapon  systems,  the  allocation 
of  space  in  ship  design  is  critical.  Modern  marine  propulsion  plants 
must  be  designed  to  provide  maximum  power  with  the  smallest, 
lightest,  as  well  as  the  most  cost  effective  equipment.  Hem 
exchangers  constitute  a  significant  portion  of  that  equipment. 

The  Naval  Postgraduate  School  (NPSi,  m  association  with  Lawci 


Tay:or  Research.  Center  and  with  supper’  from,  the  National  Science 
Foundation,  is  currently  conducting  research  tc  reduce  the  size  and 
w  eight  of  condenser  a.  Present  condensers  are  larger  than  actual.') 
required  because  of  tile  use  ca'  smooth  lubes  and  conservative  design 
margin?  The  heat -transfer  performin' m  of  a!:  condense  tube-  :• 
limited  by  the  thermal  resistances  present  or.  the  steam-side,  m  tr. 
tune  wall,  ana  on  the  water-side.  By  ieuucuq-  an\  of  these  resistance:.', 
periorrnance  is  increased  and  weight  and  size  are  decreased 


E.  CONDENSATION 

Condensation  occurs  m  two  different  modes:  fiimwisc  and 
dropwise.  Dropwise  condensation,  where  thousands  of  microscopic 
size  discrete  droplets  form  randomly  on  the  condensing  surface. 
Fiimwise  condensation  occurs  when  the  entire  condensing  surface  is 


covered  with  a  thin,  continuous  film  of  condensate  material. 
Dropwisc  is  the  preferred  mode  because  heat-transfer  takes  place 
through  a  much  thinner  effective  thickness  of  condensate  and  is 
therefore  more  efficient.  However,  the  necessary  conditions  for 
predictable,  continuous  dropwise  condensation  have  not  been 
practically  obtained.  Condenser  design  assumes  that  filmwise 
condensation  will  usually  occur,  so  that  a  conservative  margin  exists  in 
present  equipment. 

C.  HEAT -TRANSFERANALYSIS 

The  well-known  relationship  between  temperature  difference  and 
heat-transfer  rate  in  a  heat  exhanger  is: 

q=ua(lmtd>  an 

where 

O  =  heat-transfer  rate  (W». 

U  =  overall  heat  tiansle:  coefficient  (W/rn^K... 

A  =  surface  area  for  heat-transfer  consistent  with  U  t'm-n,  and 

LMTD  =  log  mean  temperature  difference  (K) 

The  overall  heat-transfer  coefficient  (U)  is  a  function  of  the  thermal 
resistances  of  the  water  side,  the  vapor  side  and  the  resistance  of  the 
tube  wall.  After  tube  material  is  selected,  the  tube  wall  resistance  is  a 
constant.  Generally,  the  thermal  resistance  on  the  water  side 
dominates.  Water-side  enhancement  techniques  include  turbulence 


promoters  and  twisted-tape  inserts,  however  these  may  come  with  an 
undesirable  increase  in  required  pumping  power  caused  by  the 
increased  internal  pressure  drop.  In  addition,  these  tube-side 
enhancement  techniques  may  be  susceptible  to  additional  fouling. 
Although  the  vapor-side  resistance  is  smaller  than  the  water-s’de 
resistance,  a  measurable  size  and  weight  benefit  may  still  result  by 
using  vapor-side  enhancement.  Vapor-side  enhancement  techniques 
include  the  use  of  low  integral  fins,  roped  tubes,  fluted  tubes,  attached 
drainage  strips,  and  dropwise  promoting  coatings.  These  methods 
have  been  reviewed  by  Marto  |Ref.  1], 

The  overall  heat-transfer  coefficient  is  commonly  based  on 

outside  surface  area  (Ao)  and  is  denoted  as  Uo,  Therefore,  upon 
neglecting  fouling  and  noncondensable  gas  resistances,  the  overall 
thermal  resistance  may  be  written  : 

l/(UoAo)  =  l/(hiAi)  +  Rw+  l/(hoAo)  (1.2) 

where 

hi  =  inside  heat-transfer  coefficient  (W/m^K), 

ho  =  outside  heat  transfer  coefficient  (W/m^-K), 

Ao  =  outside  surface  area  (m*), 

Ao»ttDrL  (1.3) 

Dt  =  fin  root  diameter  (m). 


3 


L  =  active  condensing  length  of  tube  (m),  and 
Rw  =  wall  thermal  resistance  based  on  outside  surface  area  (K/W), 
where 

Rw  =  In  l  Dr  /  Di  ]  /  (  2rtkvL)  ( 1.4) 

Di  =  inside  diameter  (m),  and 

kv  =  thermal  conductivity  of  tube  wall  material  (W/m-K). 

The  value  of  the  LMTD  in  a  condenser  depends  on  the 
characteristics  of  the  vapor/condensate  and  the  cooling  water  flow 
through  the  tube.  The  saturation  temperature  of  the  vapor  and  the 
temperature  of  the  cooling  water  source  effectively  limit  the  value: 

LMTD  =  (72  -  Ti)  /  In  [0»t  -  T i)/(T  st?  -  T2)]  (1.5) 

where 

T:  =  cooling  water  inlet  temperature  (°C), 

T 2  =  cooling  water  outlet  temperature  (°C),  and 
T si-  =  vapor  saturation  temperature  (°C). 

D.  CURRENT  TREND  S 

Modem  steam  condensers  are  equipped  with  smooth  tubes.  To 
increase  their  heat-transfer  rate,  the  addition  of  integral  fins  has  been 
investigated.  Wanniarachchi  et  al.  {Ref.  2)  reported  that  externally 
finned  tubes  have  been  widely  used  in  the  refrigeration  industry  for 


many  years,  but  that  such  tubes  were  thought  to  be  unsuitable  for  use 
in  steam  condensers.  This  belief  was  based,  in  part,  on  the 
observation  that  the  high  surface  tension  of  water  causes  liquid 
retention  (flooding)  of  the  bottom  portion  of  these  tubes.  However, 
after  recent  research  was  conducted,  their  data  and  the  data  of  other 
investigations  have  shown  that  extended  surfaces  can  provide 
significant  improvements  in  steam-side  heat-transfer.  Yau  et  al.  [Ref. 
3j  ancl  Wanniarachch:  et  al.  |Refs.  2,4j  report  that  for  condensation  of 
steam,  under  most  conditions,  the  heat-transfer  performance  of 
finned  tubes  exceeds  the  performance  expected  on  the  basis  of  area 
increase  alone. 

Marto  [Ref  5j  provides  a  comprehensive  review  of  both 
experimental  and  theoretical  developments  in  the  study  of 
condensation  on  horizontal  integral- fin  tubes.  He  states  that  no 
generalized  design  method  exists  to  accurately  predict  the  thermal 
performance  of  surface  condensers  operating  with  a  bundle  c 
horizontal  finned  tubes.  A  systematic  study  to  determine  the  effects 
of  varying  fin  geometry  and  vapor/ liquid  environment  began  a;  NP5  in 
June  1982.  Van  Petten  [Ref.  6)  provides  a  synopsis  of  that  research. 
The  majority  of  this  previous  work  has  concentrated  its  focus  on  the 
performance  and  characteristics  of  single  horizontal  tubes.  It  is  well 
known  that  condensation  in  an  actual  multi-tube  condenser  is  much 
different  than  that  which  occurs  on  a  single  tube.  Condensate 


inundation  is  one  reason  for  the  difference.  The  amount  of  published 
material  on  the  effect  of  condensate  inundation  upon  finned  tube 
performance  is  very  small  compared  to  that  of  single  tubes,  and  most 
of  it  pertains  to  low  surface  tension  fluids  used  in  refrigerant  cycles. 

E.  OBJECTIVES 

The  main  objectives  of  this  thesis  were: 

1.  Collect  steam  condensation  data  on  single  horizontal  tubes  of 
differing  diameters  for  comparison  with  previously  obtained  NPS 
results. 

2.  Modify  a  tube  being  used  in  a  concurrent  research  project  at 
Queen  Mary'  College  (University  of  London)  and  collect  steam 
condensation  data  in  an  attempt  to  resolve  an  inconsistent  trend 
previously  observed  when  comparing  enhancement  ratios  with  this 
tube  for  condensation  of  R- 113  and  steam. 

3.  Design  and  construct  modifications  to  existing  apparatus  to 
support  a  study  of  condensate  inundation  effects. 
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II.  LITERATURE  SURVEY 


A.  INTRODUCTION 

During  filmwise  condensation,  the  condensate  wets  the  solid 
surface  and  spreads  out  ia  all  directions  to  form  a  continuous  film. 
This  layer  forms  a  relatively  high  conduction  resistance  to  heat 
transfer,  as  heat  must  be  conducted  through  a  liquid  with  a  relatively 
low  thermal  conductivity.  Additionally,  this  resistance  increases  as  the 
him  tmckness  increases;  on  smooth  horizontal  tubes  him  thickness 
increases  circumferentially  from  the  top  to  the  bottom  of  the  tube.  In 
order  to  enhance  heat-transfer,  this  condensate  film  must  be  reduced. 
Thinning  of  the  condensate  film  may  be  accomplished  by  the  use  of 
tinned  surfaces,  wire-wrapped  or  roped  rubes,  and  porous  drainage 
strips.  The  focus  of  this,  and  recent  NPS,  research  has  beer,  the 
systematic  study  of  enhancement  by  the  use  of  finned  surfaces. 

During  film  condensation  on  finned  tubes,  there  are  a  large 
number  of  variables  which  can  affect  the  heat-transfer  process.  The 
three-dimensi ona!  condensate  flow'  pattern  around  the  fins  anc  oft  the 
tube,  vapor  shear  forces,  surface  tension  forces,  wall  conductor, 
effects,  and  condensate  retention  between  fins  can  all  be  important 
parameters.  The  surface  of  a  firmed  tube  during  condensation  has  two 
distinct  regions:  flooded  and  unflooded.  The  presence  of  fins 
increases  the  condensing  surface  area  and  in  the  unfiooded  region, 


surface  tension  forces  thin  the  condensate  film  over  the  fin  tips  as 
well  as  on  the  fin  flanks.  The  combined  effects  of  surface  tension  and 
gravitational  forces  then  cause  condensate  to  drain  down  into  the  fin 
root  channel,  thereby  exposing  the  maximum  amount  of  tube  surface 
for  heat-transfer.  However,  in  the  lower,  flooded  portion  of  the  tube, 
surface  tension  forces  cause  the  condensate  to  be  retained  between 
the  fins,  decreasing  the  effective  heat-transfer  area. 


The  use  of  low  integral-fin  tubes  for  enhanced  hear  transfer  with 
high  surface  tension  fluids  like  water  was  previously  considered 
impractical.  Thick  films  of  retained  condensate  on  the  lower  portion 
of  lubes  were  thought  to  effectively  eliminate  heat-transfer,  but 
Wanntarachchi.  et.  a’..  |Rel.  7]  and  Yau,  el.  al.  [Ref.  8j  have  reported 
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fully  flooded  tubes 


Apparently. 


the  flooded  region,  heat-transfer  at  the  tips  of  the  fins  is  still  effective. 
There i ore.  a  niort  complete-  understanding  of  the  processes  which 
control  heat-transfer  in  the  flooded  region  must  be  obtained. 

B.  CONDENSATE  RETENTION  ANGLE 


The  condensate  retention  angle  (¥)  has  been  defineci  as  the  angle 
measured  from  the  bottom  of  a  horizontal  finned  tube  to  the  place 
around  the  tube  circumference  where  the  condensate  film  just  fills 

the  interfin  space.  The  term,  flooding  angle  (Of)  often  appears  in  the 


literature,  ana  is  me  angle  measured  from  the  top  ot  y  horizontal 


finned  tube  to  the  place  around  the  tube  circumference  where  the 
condensate  film  just  fills  the  interfin  space.  Thus, 

'Tf=180-t  (2.1) 

A  schematic  representation  of  retention  angle  is  presented  in  Figure 

2.1. 

Van  Petten  |Ref.  6j  provides  a  detailed  account  of  retention  angle 
theory,  so  only  a  short  review  is  presented  here.  The  first 
investigation  of  retention  angle  was  conducted  by  Kata,  el.  al.  [Ref.  9,. 
They  developed  an  equation  based  on  observations  and  measurements 
taken  on  various  fluids  under  static  conditions  (no  condensation),  and 
found  that  fin  geometry  and  fluid  properties  were  controlling  factors. 
Equation  2.2  shows  the  relationship  they  derived. 

4'/ sin4-'--[o/(pfg)][(4D  i-2Dr-2s)l/[(  n/4)(D  f-Dx^Vi  (180/980, 

2  O 

v\  hen 

S'  =  condensate  retention  angle  (degrees,. 
c,  =  surface  tension  (N/np, 
g  =  gravitational  constant  (9.799  m/s£), 


pf  =  density  of  condensate  film  (kg/ms), 


Df  =  fm  tip  diameter  (m),  and 


Figure  2. 1  Schematic  of  Condensate  Retention  Angle 

on  Finned  Tube 
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s  =  interim  spacing  for  rectangular  fin  (m). 

In  the  above  expression,  as  the  ratio  of  surface  tension  to  condensate 
density  is  increased  and  as  fin  spacing  is  decreased,  retention  angle 
increases,  and  as  fin  height  and  spacing  are  kept  constant,  a  decrease 
in  root  diameter  will  produce  an  increased  condensate  retention 
angle. 

In  1981,  Rudy  and  Webb  (Ref.  10]  made  measurements  for  several 
fluids  under  both  static  and  dynamic  conditions.  They  continued  the 
Katz,  et.  a!,  relationship  and  also  concluded  that  the  condensate 
retention  angle  was  not  altered  significantly  under  dynamic  conditions 
li. e.  with  condensation).  To  account  lor  variations  in  fin  geometry. 
Rudy  and  Webb  [Ref.  1  1]  developed  the  following  relationship: 

T  =  cos*i{  1-|(2  o(2e-t))/:;D  fpfgesg  (2.8; 

whvT  v 

e  =  tin  height  (m)  and 
t  -  fin  thickness  for  rectangular  fin  (ml. 

This  expression,  shows  similar  effects  of  surface  tensio.,  tc 
condensate  density'  and  root  diameter  as  Equation  (2.2). 

A  photographic  study  was  undertaken  by  Honda,  et.  al.  [Ref.  12], 
which  confirmed  the  conclusion  of  Rudy  and  Webb  that  static  and 
dynamic  retention  angles  were  approximately  the  same.  They 
developed  the  following  relationship: 


¥  =  cos'1!  1-4  ocosJVDfp/gs] 


(2.4} 


where 

p  =  fin  rip  half  angle  (degrees). 

Rudy  and  Webb  [Ref.  13]  modified  their  previous  relationship  to 
account  for  the  general  geometry'  of  a  trapezoidal  fin: 


=  co?- 


■l{2-[(2  0(tt+2e-t  b))/t D  fplgu  te-Sbe-et  b/]i 


Where 


u  =  fin  thickness  at  tip  (m;, 
tb  =  fin  thickness  at  base  (m).  and 


sb  =  interim  spacing  at  base  (m). 

For  a  rectangular  fin,  Equation  (2.5)  reduces  to: 


t  =  cos  1  [  1-  4o/Dfpigs:  r 

which  reduces  to  Equation  (2.3)  when  (2e  -  t)  ~  2e  (i.e.,  fm  thickness 
is  small  compared  to  fir,  height),  and  identical  to  the  Honda 


relationship,  Equation  (2.4),  when  p  =  0  degrees.  As  before,  both 
Equations  (2.5)  and  (2.6)  yield  a  decreas.ng  retention  angle  as  root 
diameter  is  increased. 

Masuda  and  Rose  [Ref.  14]  conducted  a  careful  photographic  study 
of  the  condensate  flooding  process  which  revealed  that  condensate 
was  also  retained  at  the  base  of  the  fin  in  the  upper  portion  ot  the 


tube,  in  what  was  previously  termed,  the  unflooded  region.  They 
defined  several  separate  conditions  based  on  the  meniscus  profile  in 
the  interfin  spacing.  As  a  result,  Masuda  and  Rose  IRef.  141  gave  four 
different  definitions  of  flooding,  and  stated  that  Equation  (2.6)  is  valid 
for  the  most  common  definition.  Figure  2.2  is  a  schematic  of  the  four 
conditions  which  they  identified.  Nlarto  et.  al.  [Ref.  15]  observed 
similar  condensate  profiles  during  condensation  of  steam  on  wire- 
wrapped  tubes. 

I  * 

C.  HEAT -TRANSFER  PERFORMANCE  OF  A  SINGLE  TUBE 
1.  Smooth  Tube 

The  basis  tor  analysis  ot  lamina-  him  condensation  on  a 
horizontal  tube  was  formulated  by  Nusselt  [Ref.  16]  in  1016.  His 
analysis  yielded  the  well-known  relationship  for  the  heat-trenste; 

Cocliic.t : 

hNv=0./28<-K  ,D  oUf(i  s*t-T  will 

w  n  e  i  s 

k {  =  thermal  conductivity  of  the  condensate  film  (W/m-K), 
hfg  =  latent  heat  of  vaporization  (J/kg), 

p- j  =  dynamic  viscosity  of  the  condensate  film  (N-s/m2),  and 


Tv  =  wall  temperature  (°C)  (assumed  constant  around  the  tube,. 


As  Nobbs  [Ref.  17]  states,  the  Nusselt  relationship  is  subject  to  the 
following  restrictions: 

1.  the  wall  temperature  is  constant 

2.  the  flow  is  laminar  in  the  condensate  film 

3.  the  film  thickness  is  small  compared  to  normal  tube  diameters 

A.  the  fluid  properties  are  constant  within  the  condensate  film 

5.  the  heat-transfer  in  the  film  is  by  conduction 

6.  the  forces  due  to  hydrostatic  pressure,  surface  tension,  inertia 
and  vapor/liquid  interfacial  shear  are  negligible  when  compared  to 
viscous  and  gravitational  forces 

7.  the  vapor/liquid  interface  and  the  surrounding  condensate  are  at 
the  saturation  temperature 

Hopkins  [Ref.  18|,  who  conducted  research  at  NPS  on  the 
effect  of  vapor  velocity’,  details  the  developments  of  Shekriladze  and 
Gomelauri  [Ref.  19]  and  Fujii  ef.  al.  [Refs.  20  and  21]  in  modifying  the 
Nusselt  relationship  to  account  for  vapor  shear  effects.  Fujii  et.  al. 
[Ref.  201  theoretically  investigated  laminar  filmwise  condensation  of 
downward-flowing  vapor  on  a  single  horizontal  cylinder  and  they 
developed  the  following  expression  for  the  average  heat-transfer 
coefficient: 


NuRes*  =X(  1  +0.276FX  -4) 0,25  (2.8) 


where 


Nu  =  Nusseit  number  given  by  hoD<Vkf 
Re2*  =  two-phase  Reynolds  number  given  by, 


Re24=pfU«Do/H£  (2.9) 

U«»=free  stream  vapor  velocity  (m/s), 

X  =  0.9(1+ 1/G)1/3  (2.10) 

and 

G=((  VT  fkr)/(  pfhfg))((  p4tf)/(pvHv))05, 

where 

'  T'  .  nr 

Lsi  1  =  X  SAt-  i  V» 

pv  =  density  of  vapor  (kg/nh’y 

Mv  =  dynamic  viscosity  of  the  vapor  (N-s/ric;.  and 

F  =  gDcyfhfg  /(U^kfYT 
Ahe;  considering  the  more  general  case  ol  approach  velocities  in  a 
tube  bundle.  Fuju  et.  al.  [Ref.  21!  arrived  at  the  tollowmg  empirical 
correlation: 


(2. 


(2.12) 


NuRe2*  =  0.96F1/'  (for  0.03  <  F  <  600;  (2.1  A, 


2.  Finned  Tube 


Marro  [Kef.  o',  and  Van  Petten  [Ref.  6]  both  provide  an 
extensive  review  of  the  development  of  horizontal  finned  tube  heat- 
transfer,  so  just  a  brief  description  is  provided  here. 

Beatty  and  Katz  [Ref.  22]  were  the  first  to  develop  a  simple 
analytical  model  to  predict  the  condensing  heat-transfer  coefficient  on 
a  horizontal  finned  tube.  They  combined  the  Nusselt  models  for 
vertical  plates  and  smooth  horizontal  tubes  and  derived  the  following: 

hEF =0.689  ,rU.26|(Ar'Ac)(  1/D  1.3  iHfAtfAcXI/Lt0*2®)! 

(2. 1  h 

where 

F-  ( ( k  fg h  tj)/  ( (ifdT  v  i) ;  ( 2. 1 G ) 

and 


ant. 


-  nj(D 


-  Dr *>/(4D  fV 


rj f  =  fin  efficiency. 

AT  vs  =  temperature  drop  across  condensate  film  (K), 
At  =  tube  surface  area  between  fins(m^). 

Af  =  surface  area  ol  fir.  (mn  (including  Hanks  and  tip 


n  t 


Note  that  Equation  (2.15)  yields  a  heat-transfer  coefficient  which 
decreases  as  tube  diameter  is  increased. 

Rudy  and  Webb  [Ref.  10]  found  that  the  Beatty-Katz  model 
over-predicted  results  for  high  surface  tension  fluids.  They  tried  to  ust 
the  Beatty-Katz  expression  for  the  unflooded  region  only  and 
neglected  all  heat-transfer  in  the  flooded  region: 


where 


hBK*  =  (( 180-4' •/lf')hBr 


(2.19; 


hEr*  =  modified  Beatty-Katz  heat-transler  coefficient  'W/m^K; 


Owen  el.  i-\  (Kef.  23  j  concluded  that  the  Ruriy-Wehn  moat) 
underprechcted  experimental  results  because  of  the  assumption  that 
no  heat-transfer  took  place  in  the  flooded  region.  His  model 
attempted  to  include  hear  transfer  in  the  flooded  region  but  provided 
only  slightly  improved  results.  Th  :se  analyses  neglected  completely 
the  effect  that  surface  tension  forces  play  in  thinning  the  condensate 
film. 


Marto  [Ref.  6]  points  out  that  the  first  analysis  to  recognize  the 
contribution  o’  surface  tension  forces  was  Karkhu  and  Borovkov  [Rei 


24).  Webb  et.  ai.  [Ref.  25)  re-examined  the  Beatty-Katz  model  and 
concluded  that  in  the  unflooded  portion  of  the  tube,  surface  tension, 
not  gravity,  was  the  controlling  factor  in  the  movement  of  condensate. 
They  proposed  a  simplified  surface-tension-driven  model  on  the  fin 
flanks,  but  this  theoretical  result  did  not  agree  well  with  experimental 
results.  Therefore,  Webb  et  al.  [Ref.  25]  made  further  modifications 
using  the  mathematical  development  of  Adamek  (Ref.  26]  and  obtained 
results  which  were  slightly  better  than  in  their  previous  model. 

Honda  et.  al.  [Ref.  27]  have  also  recently  produced  a 
prediction  model  which  shows  considerable  promise.  However,  it 
entails  an  intensive  computer-solved  iteration  procedure  that  is 
rather  cumbersome  to  use.  In  addition,  the  probability  that  actual 
processes  occurring  on  a  flooded  tube  are  more  complex  then  can  be 
explained  in  a  one-  or  two-dimensional  conduction  model  is  apparent 
A  mere  practical,  simplifed  model,  that  is  adaptable  tc  universal 
design  use.  is  required. 

D.  HEAT  TRANSFER  PERFORMANCE  INCLUDING 

CONDENSATE  INUNDATION 

Morrison  [Ref.  28]  and  Brower  [Ref.  29]  were  involved  in  previous 
NPS  research  on  the  effects  of  condensate  inundation  on  enhanced 
condenser  tubing.  They  provide  a  comprehensive  review  of  the 
pertinent  literature  up  through  1985. 


1.  Smooth  Tube 


Jakob  [Ref,  30]  extended  the  Nusselt  analysis  to  film 
condensation  on  a  vertical  in-line  column  of  horizontal  tubes  by 
assuming  that  all  the  condensate  from  a  tube  drains  as  a  laminar  sheet 
onto  the  tube  beneath.  For  this  ideal  condition,  the  average  coefficient 
for  a  vertical  column  of  N  tubes  was  predicted  as: 

hw  =  0.728[(kf3pfahfgg)/(D  oNjtf(AT  f))j  0-25  (2.20) 

When  Equations  (2.7)  and  (2.20)  are  combined,  the  Nusselt  theory  for 
the  average  coefficient  for  a  column  of  N  tubes  compared  to  the 

coefficient  of  the  top  tube  (hi)  gives: 

hu/hi  =  N*1/4  (2.21) 

Kern  [Ref.  31)  realized  that  condensate  does  not  flow  in  a 
laminar  sheet,  but  by  discrete  droplets  or  columns.  Ke  proposed  a 
less  conservative  expression  to  account  for  the  disturbances  in  the 
condensate  film.  The  Kern  relationship  is: 

hN/hi  =  N-1/6  (2.22) 

Eissenberg  [I  ef.  32]  conducted  experimentation  with  a 
staggered  tube  bundle  which  included  the  effects  of  steam  velocity, 
condensate  inundation,  and  noncondensable  gases  on  the 
condensation  heat-transfer  coefficient.  As  a  result  of  his 
investigations,  he  developed  a  side-drainage  model.  In  this  model. 


condensate  may  drain  to  tubes  that  are  not  directly  below  one  another. 
This  side  drainage  could  be  influenced  by 

1.  bundle  orientation  -  triangular  spaced  arrays  are  more  susceptible 
to  side  drainage  than  are  inline  tube  arrays 

2.  tube  pitch  -  as  it  becomes  smaller,  side  drainage  should  occur 
more  frequently 

3.  momentum  -  the  greater  the  horizontal  momentum  of  a  drop 
leaving  the  tube,  the  greater  is  its  tendency  to  undergo  side 
drainage 

4.  vapor  velocity  and  direction  -  the  more  horizontal  it  is,  the  more 
lateral  momentum  it  will  impart 

5.  tube  alignment  -  a  tube  misaligned  in  a  bundle  may  receive  and 
then  release  a  greater  or  lesser  amount  of  condensate 

If  the  condensate  drains  purely  by  side  drainage,  Eissenberg's  model 
predicted: 

hN/hNu  =  0.60+0.42N"1/4  (2.24) 

Figures  Z2  a,b,c  depict  the  three  models  theorized  by 
Nusselt,  Kern,  and  Eissenberg,  respectively.  Figure  2.3  is  a  graphical 
representation  of  the  three  theories  with  a  superimposed  spread  of 
experimental  data. 

Nobbs  (Ref.  17),  although  unable  to  find  a  general  correlation 
for  his  experimentation  on  the  combined  effects  of  vapor  velocity  and 


condensate  inundation  in  a  tube  bundle,  did  supply  the  following 
conclusions: 

•  vapor  velocity  increases  the  condensate  heat-transfer  coefficient 
on  both  inundated  and  umndated  tubes 

•  the  effect  of  inundation  generally  reduces  the  heat-transfer 
coefficient 

•  the  condensate  drainage  path  is  often  not  vertically  downwards, 

but  jv  m  a  diagonal  direction,  particularly  in  triangular-shaped 
arrays 

Because  experimental  results  vary  considerably,  results  have 
been  fitted  to  an  equation  using  the  following  form: 

hN/hi  =  N-S  (2.25) 

where  S  is  determined  empirically.  Reported  values  oi  S  vary  in  the 
range  or  0.07  to  0.20. 

2.  Finned  Tubes 

Very  little  data  has  been  published  for  finned  tubes,  most  of 
which  details  refrigerants  and  small  tube  bundles.  Marto  [Ref.  5? 
provides  an  accounting  of  published  works  which  indicate  that 
inundation  effects  are  less  for  finned  tubes.  Most  recently.  Honda 
et,  el  jRef.  33)  performed  experiments  using  R-113,  methanol,  and 
normal  propanol  to  study  the  characteristics  of  falling  condensate  on  a 
vertical  column  of  horizontal  low  finned  tubes.  They  reported  the 
observation  of  four  major  flow  modes:  droplet,  column,  column  and 
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Figure  2.3  The  Effect  of  Condensate  Inundation  on 
Smooth  Tube  Performance 
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sheer,  and  sheer  as  inundarion  flow  rare  was  increased.  Flow  mode 
transition  was  associated  with  the  value  of  a  dimensionless  group,  K: 

K  =  r(g/pi)°'25/a0'75  (2.26) 

where 

T  =  liquid  film  flow  rate  per  unit  length  from  one  side  of  the 
tube  (kg/m-s). 

Photographs,  as  well  as  visual  observations,  showed  that  at  relatively 
low  rates  of  inundation,  liquid  droplets  or  columns  fell  from  an  upper 
tube  onto  a  lower  one  and  flowed  around  the  tube  circumference 
through  a  number  of  grooves  between  adjacent  fins.  Thus,  the  rest  of 
the  tube  surface  was  not  affected  by  the  falling  liquid.  Therfore,  at 
these  low  inundation  rates,  heat-transfer  performance  of  finned  tubes 
w'ould  be  expected  to  be  less  affected  than  the  performance  of 
smooth  tubes. 

Honda  et.  al.  [Ref.  34]  then  made  heat-transfer  measurements 
for  the  film  condensation  of  R- 113  on  in-line  bundles  of  both  flat-side- 
finned  (rectangular- finned)  tubes  and  three-dimensional  finned  tubes 
(i.e.  spines).  The  effect  of  vapor  shear  and  condensate  inundation 
were  both  examined.  They  found  that  the  flooding  level  of  condensate 
for  the  unaffected  region  of  the  flat-sided-finned  tubes  was 
independent  of  vapor  velocity.  Their  data  demonstrated  that  heat- 
transfer  enhancement  due  to  vapor  shear  was  considerably  smaller  for 


finned  tubes  than  for  smooth  tubes,  and  that  condensate  inundation 
was  more  significant  for  the  three-dimensional-fin  tube  than  for  the 
flat-sided-finned  tube. 


III.  APPARATUS  AND  TUBES  TESTED 


A.  DESCRIPTION  OF  APPARATUS 


The  test  apparatus  used  to  conduct  the  single  tube  studies  was 
identical  to  that  utilized  by  Van  Petten  [Ref.  5j,  and  is  shown  in  Figure 
3.1.  It  consisted  of  a  boiler,  glass  piping,  auxiliary  condenser,  purging 
system,  and  test  section.  The  boiler  was  made  from  a  0.3048  m 


diameter  Pvrex  glass  section  circumferentially  fitted  from  below  with 
ten  4k\V.  440V,  Watlou  immersion  hearers.  Vapor  produced  there 
flowed  upward  through  a  2.44  m  long  section  of  0.152  rn  Pyrcx  glass 
piping;  it  negotiated  a  180  degree  turn  and  then  travelled  downward 
through  a  1.52  m  long  section  of  Pvrex  glass  piping  before  arriving  a; 
the  test  section.  The  test  section,  Figure  3.2.  was  fabricated  from 
stainless  steel  with  nvlon  and  Teflon  fitting*-  to  support  the  tes1  tube 


and  provide  cooling  water  connections.  Cooling  water  to  the  test  tube 
was  pumped  from  a  sump  tank  and  its  flow  rate  w'as  controlled  by  a 
t hr oi tie  vaive  at  the  flowmeter  inlet.  A  continuous  flow  c;  filtered 
tapwater  w>as  supplied  to  the  water  sump.  A  manometer  and 
thermocouple  placed  just  above  the  test  tube  provided  vapor  pressure 
and  temperature  readings.  A  circular  viewing  port  allowed  visual 
observation  of  the  condensing  process.  The  auxiliary'  condenser, 
located  directly  below  the  test  section,  consisted  of  two  concentrically 
wound  copper  coils  inside  Pvrex  glass  piping.  Auxiliary  cooling  water 
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DRAIN  FILL  /  DRAIN  VALVE 


Figure  3.2  Schematic  of  Test  Section  ( insert  Removed  ) 


29 


was  supplied  by  tap  water  thrcttiea  at  the  inlet  of  a  flowmeter.  The 
purging  system  consisted  cf  copper  tubing  ,  a  copper  coil.  Plexiglass 
cylinder,  and  a  vacuum  pump.  Gases  drawn  by  the  pump  left  the 
auxiliary  condenser  through  the  copper  tubing  and  flowed  into  the 
copper  coil  which  was  situated  inside  the  cooling  water  sump.  It  was 
hoped  that  heat  exchange  here  would  condense  any  residual  working 
vapor  before  it  reached  the  operating  vacuum  pump.  The 
noncondensable  gases  and  residual  working  fluid  were  collected  in  the 
Plexiglass  cylinder  and  the  fluid  was  emptied  using  a  drain  valve  a: 
the  conclusion  of  each  data  set.  The  gases  were  sucked  to  the  top  ol 
the  container  and  draum  into  the  vacuum  pump.  Figure  3.3  is  a 
schematic  of  the  purging  system.  Other  system  details  can  ne  found  m 
(F P  ]J-<  no  np  07  no; 
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B.  MODIFICATIONS  TO  APPARATUS  FOR  INUNDATION  TESTS 

Modifications  were  made  to  the  apparatus  to  allow  to-  the  s’udy  r  : 
condensate  inundation.  The  test  section  used  by  Hopkins  [Ref.  R 
was  used  after  several  modifications  were  installed.  This  test  section 
was  designed  to  support  the  installation  and  visual  obervation  oi 
several  in-line  tubes.  Removable  tube  sheets  manufactured  trom 
Teflon  (with  aluminum  backing  plates  for  rigid  fastening  support  and 
scaling)  allowed  for  the  possible  installation  of  up  to  five  tubes. 
Hopkins  conducted  his  investigation  using  a  single  tube.  He  installed  a 


Figure  3.3  Schematic  of  Purging  System 
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removable  converging  nozzle  arrangement  in  order  to  increase  the 
vapor  velocity  through  the  test  section.  This  section  was  removed, 
leaving  a  test  section  almost  identical  to  the  single-tube  apparatus. 
New'  inlet  and  outlet  tube  blanks  were  manufactured  to  support  a  pitch 
to  diameter  ratio  of  1.25  for  the  four  tubes  that  would  be  installed.  A 
laser  sighting  procedure  was  utilized  to  ensure  proper  vertical  and 
horizontal  alignment  of  the  four  tubes  within  the  test  section.  Figure 
3.4  is  a  schematic  of  the  test  section  (which  looks  at  rhr  tube  blank 
side).  Figure  3.5  is  a  cross-sectional  view  of  the  test  section  showing 
the  installed  tubes  .  The  bottom  tube  was  the  active  tube  upon  which 
condensation  took  place  and  for  which  cooling  water  temperature 
increase  measurements  were  taken.  The  tube  above  it  matched  the 


geometric  parameters  exactly  to  ensure  that  the  flow  mode  of  the 
falling  condensate  w'ould  not  be  altered  because  of  geometric 
configurational  differences.  The  ton  two  lubes  were  inundation  tubes 


connecter,  in  r.  condensate  return  system,  me  oripma.  oesm- 


consisted  ot  a  larger,  constant  speed  pump  throttled  through  rwe 
flowmeters,  but  operation  of  that  configuration  was  plagued  by 
cavitation  of  the  pump.  The  second  system,  designed  only  fev  flow 
visualization,  W'as  capable  of  supporting  a  single  inundation  tube.  The 
pump,  throttling  valves  and  flowmeters  were  replaced  with  copper 
and  tvgon  tubing,  and  a  condensate  pump  with  a  controller  that 
accurately  varied  its  output  through  a  calibrated  range.  Condensate 


3.4  Schematic  of  Test  Section  Showing  Tube  Sheet 
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Figure  3.5  Cross  -  Sectional  View  of  Test  Section 
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was  suctioned  from  the  condensate  drain  line  by  the  condensate 
pump  and  sent  upward  toward  the  test  section.  There  it  was  split  and 
approached  the  bottom  inundation  tube  from  both  ends  (on  an  equal 
pressure  drop  basis).  The  inundation  tubes,  pictured  schematically  in 
Figure  3.6,  were  manufactured  so  that  their  internal  and  external 
diameters  matched  the  finned  tubes'  internal  and  root  diameters, 
respectively.  Condensate  was  forced  outward  through  a  series  of  27 
1.0  mm  diameter  holes,  drilled  on  5  mm  centers,  in  the  top  of  the 
tube.  An  umbrella/ cap  mounted  longitudinally  atop  of  the  tube 
deflected  all  the  condensate  downward  along  the  entire  tube's  surface 
where  it  could  then  freely  flow  down  to  the  the  tubes  below.  This  was 
needed  to  prevent  geysering  of  the  liquid  at  high  inundation  flow 
rates.  Figure  3.7  is  a  photograph  of  both  inundation  tubes.  The  top 
tube  is  viewed  from  the  side,  clearly  showing  the  umbrella  in  place. 
The  lower  tube  is  viewed  from  the  top.  A  schematic  diagram  of  the 
modified  apparatus,  with  the  original  condensate  return  system, 
appears  as  Figure  3.8. 

C.  TUBES  TESTED 

Two  smooth  tubes  and  14  finned  tubes  were  to  be  tested  with 
steam  as  the  working  fluid.  All  tubes  were  manufactured  from  copper. 
Of  the  finned  tubes,  12  of  them  consisted  of  two  families  of  six  tubes 
each  whose  root  diameters  differed.  Those  with  a  root  diameter  of 
12.7  mm  were  termed  small  tubes,  and  those  with  a  root  diameter  of 
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Figure  3.8  Schematic  of  Modified  Test  Apparatus 
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19.05  mm  were  termed  medium  tubes.  Van  Petten  [Ref.  5]  tested  an 
additional  family  of  tubes  whose  root  diameter  was  25.0  mm.  These 
were  known  as  the  large  tubes.  Each  tube  in  a  family  matched  its 
corresponding  member  in  the  other  family,  fin  height,  thickness,  and 
spacing  were  identical  These  finned  tubes  had  a  fin  height  of  1  mm 
and  a  fin  thickness  of  1  mm,  with  fin  spacings  of  0.25.  0.50,  1.0,  1.5, 
2.0  and  4.0  mm.  Each  smooth  tube's  outside  diameter  matched  the 
root  diameter  of  its  corresponding  family.  Figure  3.9  is  a  photograph 
of  the  medium-diameter  tube  family  and  Figure  3.10  is  a  photograph 
of  small-  and  medium-diameter  tubes  with  a  1.50  mm  fin  spacing. 

The  remaining  two  finned  tubes  were  the  basis  for  a  comparison  of 
results  beng  obtained  in  a  collaborative  experimental  program  with 
Queen  Mary  College  (QMC)  of  the  University  of  London.  Van  Petten 
[Ref.  6]  conducted  preliminary  research  using  a  tube  which  had  been 
manufactured  to  meet  the  geometric  specifications  of  the  a  tube  beinf, 
used  in  the  QMC  study  (termed  the  QMCNPS  tube).  However,  as  he 
reported,  data  taken  for  R-113  and  steam  appeared  to  6how 
unexplainable  inconsistencies.  Therefore,  a  tube  that  was  used  in  the 
QMC  study  of  Masuda  and  Rose  [Ref.  40)  was  obtained  from  QMC  and 
was  then  modified  for  use  in  the  NPS  apparatus  (termed  the  QMC 
tube).  Tube  parameters  are  listed  in  Table  3:1.  As  denoted  in  Tabic 
3.1,  and  elaborated  upon  in  Table  3.2  there  were  some  small 
differences  between  the  QMC  and  QMCNPS  tubes.  Because  of  it s 
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DIMENSIONS  OF  TUBES  TESTED 


Tube 

Dr  (mm) 

Di(mm) 

s  (mm) 

t  (mm) 

e  (mm) 

Small  Smooth 

12.70 

9.53 

- 

- 

- 

Small  Finned 

2 

12.70 

9.53 

0.25 

1.0 

1.0 

3 

12.70 

9.53 

0.50 

1.0 

i.O 

4 

12,70 

9.53 

1.0 

1.0 

1.0 

5 

1 2.70 

9.53 

1.5 

1.0 

1.0 

6 

12.7o 

9.53 

2.0 

1.0 

1.0 

7 

12.70 

9.53 

4.0 

1.0 

*  .( ' 

QMC' 

12.60 

9.84 

1.0 

0.5 

1.59 

QMCNPS’ 

12.70 

9.53 

1.0 

o.r> 

1.5:’ 

Medium  Smooth 

19.05 

12.70 

- 

- 

- 

Medium  Finned 

8 

19.05 

12.70 

0.25 

1.0 

k  .  V  ' 

o 

19.05 

12.70 

0.50 

1.0 

t  n 

A  •  8* 

JO 

19.04 

}  2.70 

l.(i 

* 

!  • 

X  X 

19.05 

12.70 

1.5 

1.0 

1.0 

i  n 

X  -- 

19.05 

12.70 

2.0 

1.0 

1  (  ? 

i  v 

19.05 

12.70 

4.0 

1 .!. 1 

i  .*  ? 

*  Average.-  valuers.  seeTable  3.2  for  details 

TABLE  3.2 

COMPARISON  OF  QMC  AND  QMCNPS  TUBES 


QMC  QMCNPS 


Inside  diameter  (mm) 
Outside  diameter  (mm) 
Overall  length  (mm) 
Condensing  Length  (mm) 
Fin  Thickness  (mm) 

Fin  Spacing  (mm) 

Fin  Height  (mm, 


9.84  9.53 

12.60  12T0 

167.5  228.0 

100.0  133.4 

0.38  -  0.7 1  0.48-0.5  1 

0.94  -  1.02  0.94-1.0'. 

1.59  1.59 


-non ur  overa.i  iength.  the  QMC  rube  had  ro  be  modified  in  order  ro 
be  rested  in  the  NPS  apparatus.  Adaptors  were  designed  for  each  end 
to  match  the  lengths  of  the  inlet  and  outlet  sections  of  the  standard 
NPS  tube  and  then  soft-soldered  in  place.  The  QMC  tube 
also  has  a  shorter  condensing  length  so  insulating  Teflon 
sleeves  were  placed  at  either  side  of  the  active  condensing  finned  area 
to  prevent  condensation  on  the  remaining  exposed  area.  Finally  it  also 
had  a  small  (~  1  mm  wide)  longitudinal  slot  machined  in  the  fins  to 
accept  a  porous  drainage  plate.  When  tested,  the  tube  was  oriented  so 
that  this  slot  was  at  the  bottom.  Figures  3.11  and  3.12  are 
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the  two  tubes,  before  and  after  modifications  to 


QMC  tube,  respectively.  An  enlarged  photograph  of  the  tube  surlaces, 


r  ifcwic  3. 13.  along  with  the  measurements  presented  in  Table  3.2 


clearly  shew  differences  in  the  machined  quality'  of  the  two  tubes 


D.  INSERTS 

/-\s  mentioned  previously,  the  water-side  thermal  resistance  can 
contribute  substantially  tc  the  overall  thermal  resistance.  Search  [Ref. 
41J  reported  that  as  much  as  56%  of  the  overall  resistance  could  be 
attributed  to  the  water-side  in  a  Naval  condenser.  Therefore,  in 
order  to  improve  the  accuracy  in  calculating  that  outside  coefficient, 
Georgiadis  [Ref.  351  designed  and  analyzed  the  use  of  spiral  inserts  for 
the  water-side,  shown  in  Figure  3.14. 


Figure  3 .13  Enlargement  of  Machined  Surfaces, 

QMC  (Top)  anti  QMCNP3  ( Bottom 


In  rhe  absence  of  inside  enhancement,  the  tube-side  resistance  is 
substantially  greater  than  that  of  that  vapor-side.  Therfore,  an  error  in 
the  inside  resistance  results  in  a  significantly  larger  error  on  the 
outside  resistance.  Internal  enhancement  also  acts  to  reduce 
circumferential  wall  temperature  variation  and  increases  the  heat 
transfer  rate  for  a  given  cooling  water  velocity,  which  thereby  reduces 
the  relative  error  in  total  heat-transfer.  The  spiral  insert  greatly 
reduces  thermal  entrance  effects  and  their  associated  uncertainties. 

E.  SYSTEM  INTEGRITY 

System  integrity  describes  the  vacuum-tightness  of  the  test 
apparatus.  The  system  had  been  inoperative  for  a  period  of 
approximately  six  months  before  this  investigation,  so  all  of  the 
system's  rubber  gaskets  were  replaced  before  operations  commenced. 
A  vacuum  test  was  periodically  conducted  by  using  the  vacuum-purge 
system  to  bring  the  system  pressure  down  to  about  20  mm  of  mercury. 
The  system  was  secured  and  the  vacuum  system  was  turned  off.  The 
manometer  was  closely  monitored  over  a  24  hour  period  to  ensure 
that  any  increase  in  system  pressure  was  within  prescribed  limits.  A 
leak  rate  from  12-20  mm  of  mercury  over  a  24  hour  period  was 
typical  (in  earlier  investigations,  a  leak  rate  as  low  as  2  mm  of  mercury 
was  sometimes  obtained). 

The  modification  of  the  apparatus  to  the  inundation  scheme 
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introduced  a  whole  new  series  of  potential  vacuum  leaks.  To  minimize 
these  leaks,  solder  joints  instead  of  fittings  were  used  wherever 
possible,  but  even  then  the  subsequent  process  of  eliminating  leaks 
was  long  and  laborious.  A  constant  source  of  problems  was  the  seal 
around  the  tube  sheets,  this  was  finally  solved  with  the  installation  of 
larger  O-rings  and  a  thicker  outlet  aluminum  backing  plate  whose 
surface  was  milled  flat  to  ensure  a  uniform  seal  of  the  teflon  tube 
sheer. 

The  presence  of  noncondensable  gases  during  the  collection  of 
data  was  a  definite  indicator  of  vacuum  leaks.  The  data  reduction 
software  prompted  the  operator  to  check  for  the  presence  of 
noncondensables  at  the  beginning  and  end  of  each  data  run.  A 
mercury-in-glass  manometer  measured  the  absolute  pressure  in  the 
test  section  and  a  thermocouple  measured  steam  temperature.  The 
saturation  temperature  corresponding  to  the  measured  pressure,  and 
the  Gibbs-Dalton  ideal-gas-mixture  relations  were  used  to  compute  the 
noncondensable  gas  (asumed  to  be  air)  concentration.  The  data 
reduction  software  prompted  the  operator  to  check  for  the  presence 
of  noncondensables  at  each  datum  and  issued  a  warning  if  the  the 
percentage  exceeded  programmed  limits. 
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IV.  SYSTEM  OPERATION  AND  DATA  REDUCTION 


A.  SYSTEM  PROCEDURES 

Copper  has  a  tendency  to  condense  steam  in  the  dropwise  mode. 
In  order  to  ensure  that  only  filmwise  condensation  takes  place, 
previous  investigators  devised  a  procedure  for  the  preparation  of  the 
tube  surface.  The  tube's  external  and  interna!  surfaces  are  first 
cleaned  with  a  soft  brush  and  tapwater.  The  tube  is  then  suspended 
above  a  steam  bath  and  a  coating  of  a  heated  mixture  containing  50% 
ethyl  alcohol  and  50%  sodium  hydroxide  is  applied  to  the  tube  surface. 
For  an  initial  coating,  the  mixture  is  applied  once  every  ten  minutes 
for  an  hour.  If  a  tube  has  been  previously  treated,  the  mixture  is 
applied  once  every  five  minutes  over  a  20  minute  period.  This 
procedure  deposits  an  ultra-thin,  black  oxide  layer  on  the  outside 
surface  of  the  tube.  The  tube  is  then  removed  from  the  steam  bath 
and  thoroughly  rinsed  with  distilled  water  prior  to  its  installation  in 
the  test  section.  The  oxide  film  gives  the  surface  very  high  wetting 
characteristics  and  its  thermal  resistance  is  insignificant.  System 
startup  and  general  operating  procedures  are  detailed  in  Appendix  A. 

System  pressure,  measured  using  a  mercury-in-glass  manometer, 
ana  temperature,  measured  with  the  voltage  reading  of  the  vapor 
thermocouple,  were  closely  monitored  in  order  to  maintain  the 


desired  operating  conditions.  Table  4.1  lists  the  operating  conditions 
used  during  the  investigation: 

TABLE  4.1 

OPERATING  CONDITIONS 


Pressure  Conditon 

Thermocouple 

Measured 

Measured 

Reading 

Temperature 

Pressure 

(microvolts) 

(°C) 

(mm  Hg) 

Vacuum 

1977 

48 

85 

Atmospheric 

4277 

100 

765 

Once  the  operating  temperature  was  established  (about  45  minutes 
after  start-up),  sample  data  were  taken  to  check  for  non-condensables 
and  desired  vapor  velocity'.  Adjustments  to  system  temperature  were 
made  by  manually  controlling  the  amount  of  cooling  water  through  the 
auxiliary'  condenser.  The  viewing  window  allowed  direct  observation 
of  the  condensing  process.  Cooling  W'ater  flow  rate  was  manually 
emered  for  each  data  run.  Eight  different  rates  w'ere  employed  to 
obtain  the  largest  possible  range  of  cooling  water  temperature  rise.  To 
demonstrate  repeatability',  multiple  data  runs  were  taken  and  results 
compared. 

B.  DATA  REDUCTION 

As  discussed  in  Chapter  I,  the  total  resistance  to  heat  transfer 
consists  of  the  sum  of  the  resistances  of  water-side,  wall,  and  vapor- 
side.  The  w'ater-side  (inside)  and  vapor-side  (outside)  resistances  are 


convective  in  nature  and  can  be  represented  as: 

Ri  =  l/(hiAi)  (4.1) 

and 

Ro=  l/(hcAo)  (4.2) 

where 

Ri  =  inside  resistance  to  heat  transfer  (K/W), 


and 

Ro  =  outside  resistance  to  heat  transfer  (K/W). 


The  calculation  of  Ai  must  include  not  only  the  inside  surace  area 
corresponding  to  the  active  condensing  length,  but  also  the  inlet  and 
outlet  portions  of  the  tube  which  lie  inside  the  nylon  bushings  and  act 
as  fins  to  conduct  heat  away.  Therefore, 

Ai  =  itDi(L  +  Lir)i  +  hzr\2),  (4.3) 

and 


Ao  =  TlDrL, 


(4.4) 


Ai  =  effective  inside  area  (m2), 

Lt  =  length  of  inlet  portion  of  tube  (m), 


r]i  =  fin  efficiency  for  inlet  portion  of  tube, 

L2  =  length  of  outlet  portion  of  tube  (m),  and 
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qs  =  fin  efficiency  for  outlet  portion  of  tube. 

Combining  Equations  (1.2),  (4.1),  and  (4.2)  gives  the  following 
expression; 

1/U  oAo  =  Ro  +  Ri  +  Rw  (4.5) 

To  determine  the  overall  heat-transfer  coefficient,  the  heat- 
transfer  rate  must  be  calculated.  This  was  accomplished  by  measuring 
the  inlet  and  outlet  temperatures  and  the  mass  flow  rate  of  the 
cooling  water  through  the  tube.  A  correction  factor  was  applied  to 
the  ouTiei  temperature  to  account  for  the  effect  of  viscous  heating  ol 
the  cooling  water  due  tc  mixing  and  to  the  induced  swirl  from  the 
insert  Properties  of  the  cooling  water  were  calculated  based  on  an 
average  of  the  inlet  and  outlet  temperatures.  Quartz  thermometer 
readings  were  used  to  determine  the  inlet  and  oulet  temperatures.  A 
calibrated  flowmeter  provided  the  mass  flow  rate.  An  energy  balance 
gave; 

Q  =  rhcpAT  (4.6) 

where 

m  =  mass  flow  rate  (kg/s)  and 

cp  =  specific  heat  at  constant  press  ire  (J/kg-K). 

The  velocity  of  the  cooling  water  through  the  tube  was  determined  by 
the  simple  relationship; 
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V  =  m/(pcAc) 


(4.7) 


where 

V  =  test  tube  average  cooling  water  velocity  (m/s), 


pc  =  test  tube  cooling  water  density  (kg/m3),  and 


Ac  =  cross-sectional  flow  area  of  tube  (equal  to  7tDi2/4)  (m2). 

With  the  value  of  the  heat  transfer  rate  known,  and  using 
Equations  (1.1)  and  (1.4),  the  overall  thermal  resistance  is  found. 
Since,  the  waP  resistance  can  be  calculated  by  Equation  (1.4),  the  two 
unknowns  remaining  in  Equation  (4.5)  are  the  values  for  the  inside 

and  outside  heat  transfer  coefficients,  which  appear  in  Ri  and  Ro  , 

respectively.  To  obtain  these  values,  a  Modified-Wilson-Plot  technique 
was  employed.  The  alternative  was  the  use  of  instrumented  fmned 
tubes  fined  with  wail  thermocouples.  Past  investigations  used 
instrumented  smooth  tubes  and  subtracted  the  measured  inside 


to  get  the  outside  resistance.  A  comparison  of  the 


various  methods  has  shown  that  the  Modified-Wilson-Plot  method. 


when  done  properly,  yielded  r  suits  consistent  with  those  obtained  by 
the  costlier,  more  time-consuming  method  of  direct-reading, 
instrumented  tubes. 


The  Modified-Wilson-Plot  method  requires  that  the  form  for  both 
inside  and  outside  heat  transfer  coefficients  be  known.  To  maximize 
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the  accuracy,  the  inside  and  outside  resistances  should  be  of  equal 
magnitude,  A  spiral  insert  was  used  to  reduce  the  inside  resistance. 
During  this,  and  past,  NPS  studies,  a  Sieder-Tate  relationship  was 
used  for  the  inside  heat  transfer  coefficient: 


hi  =  Ci(kc/Di)Re0,8Pr0,333(pc/pw)0,14  =  CiQ  (4.8) 

where 

Ci  =  Sieder-Tate  type  coefficient, 

kc  =  thermal  conductivity  of  cooling  water  (W/m-K), 

Re  =  Reynolds  number  of  cooling  water  (pVDi/pw), 

Pr  =  Prandtl  number  of  cooling  water, 

pc  =dynamic  viscosity  of  cooling  water,  evaluated  at  av;rage 

cooling  water  temperature  (N-s/m2),  and 

p  w  =  dynamic  viscosity  of  cooling  water,  evaluated  at  wall 

temperature  (N-s/m*). 

A  Nusselt-type  equation  was  used  to  determine  the  outside  heat- 
transfer  coefficient: 

ho  =  <x[(k  fSp^ghf g)/(  pfDoq))  1/3«  ocF*  (4.9) 

where 


a  =  dimensionless  coefficient,  and 


q  =  heat  flux  based  on  outside  area  (W/m2). 

q  =  Q/A0  (4.10) 

After  the  proper  substitutions  of  Equations  (4.8)  and  (4.9)  into 
(4,5),  the  result  becomes: 

I(l/U  o)  -  Rv]F*  =  I(AoF*)/(AiCiQ)]  +  ( 1/oc)  (4.11) 

This  is  a  linear  equation  with  two  unknowns;  Ci  and  ot.  To  solve  for 
them  we  make  the  following  substitutions: 

r=  |(l/Uo)  -  Rw]F*.  (4.12) 

|(AoO/(AiO;j.  (4.13) 

Ci  =  (1/m).  (4.14) 

and 

cc  =  ( 1/b).  (4.15) 

to  arrive  at  the  following  linear  expression: 

X=  mA  +  b  (4.16) 

A  least-squares  fit  of  Equation  (4.16)  was  then  made  to  determine  <x 
and  Ci  simultaneously.  An  iterative  procedure  was  required  to 

calculate  both  F*  and  Q,  because  of  their  temperature  dependences. 

To  calculate  enhancement  ratios,  the  following  analysis  was  used. 


56 


starting  with  this  relationship; 


Aftersubstituting 


q  =  aATfn 


q  =  ho  ATf 

the  resulting  relationship  is: 


(4.17) 


(4.18) 


ho  =  aATf^-l  (4.19) 

Nusselt  theory  predicts  a  value  of  0.75  for  n.  Therfore,  n  was  set  to 
0.75  so  that  the  enhancement  ratio,  based  on  constant  heat  flux  could 
be  determined 


(hof/hos)  =  (af/az)  =  |(ctfFf*)/(  <xsFs*)]  (4.20) 

where  subscripts  f  and  s  denote  finned  tube  and  smooth  tube, 
repectiveiy. 

hof  =outside  heat-transfer  coefficient  for  finned  tube  (W/m^K). 

and 

hos  =outside  heat-transfer  coefficient  for  smooth  tube  (W/m^-K). 

If  the  heat  flux  across  the  condensate  film  is  held  constant  for  both 

finned  and  smooth  tubes  then  their  values  of  F*  will  be  equal  reducing 
Equation  (4.20)  to: 


Sq  =  (0Cf/<Xs) 


(4.21) 


The  values  for  otf  and  <xs  were  taken  directly  from  the  Modified- 
Wilson-Plot  analysis  to  calculate  the  enhancement  ratio.  For  constant 
vapor-side  temperature  drop,  Equation  (4.21)  can  be  replaced  by: 

®  at  =  (<x$/<x  s)3/4  (4.22) 

Once  determined,  by  using  Equation  (4.14),  the  value  for  the 
Sieder-Tate  coefficient  was  used  to  calculate  the  inside  heat-transfer 
coefficient  and  the  outside  heat-transfer  then  was  calculated  as: 

ho  =  U/Kl/Uo)  -  Ao/(hiAi)  -  AoRv)]}  (4.23) 

When  graphically  presenting  the  data,  least-squares  fits  of  Equation 
(4.19),  using  values  of  ho  calculated  by  Equation  (4.23),  were  utilized. 

The  effect  of  tube  diameter  appears  only  in  F*  and  therefore  a 
should  be  independent  of  the  tube  diameter.  Previously,  Van  Petten 
[Re f.  6]  calculated  enhancements  for  both  medium  and  small  tubes 

based  on  an  a  he  obtained  for  a  medium  smooth  tube.  He  did 
recommend  that  a  small  smooth  diameter  be  manufactured  to  confirm 
the  theoretical  conclusion.  A  small  smooth  tube  was  fabricated  during 

this  thesis,  and  after  testing,  an  cx.  which  differed  slightly  from  the 
medium  tube  was  obtained.  Specific  results  are  presented  in  the 
following  chapter. 

The  additional  calculations  required  for  the  inundation  data  were 


58 


minimal.  In  order  to  calculate  the  dimensionless  parameter,  K,  with 
which  Honda  et.  al.  [Ref.  35j  predicted  flow  mode  transition,  the 
inundation  flow  rate  and  the  mass  condensation  rate  were  required. 
Inundation  flow  rate,  as  measured  with  calibrated  condensate  pump 
controller  settings,  was  manually  entered  into  the  program.  The  mass 
condensation  rate  was  evaluated  as: 

me  =  Q/hfg  (4.24) 

where 

m-  =  mass  condensate  rate  (kg/s), 
r  could  then  be  calulated: 

r  =  (mc  +  mi)/(2L)  (4.25) 

where 

m i  =  inundation  flow  rate  (kg/s). 

K  could  then  be  calculated  using  Equation  (2.26).  Visual  observations 
w'ere  made  as  inundation  rate  was  increased  in  an  attempt  to  observe 
the  predicted  flow'  modes  of  droplet,  column,  column-sheet,  and 
sheet. 

C.  PROGRAM  MODIFICATIONS 

The  modifications  to  the  program  utilized  by  Van  Petten  [Ref.  6] 
included. 

1  acuir- K>nai  tube  diameter  variable  was  added  lot  the  Cyivl C  tube 


2.  exponent  for  Reynolds  number  in  the  Sieder-Tate  equation  was 
allowed  to  vary  at  operator’s  discretion 

3.  an  alternate  relationship  for  the  inside  heat-transfer  coefficient 
was  provided  as  an  option  to  the  operator 

4.  <x  for  small  smooth  tube  was  added  and  used  to  calculate 
enhancement  ratios  for  small,  QMC,  and  QMCNPS  tubes 

5.  required  inundation  inputs  and  calculations  were  incorporated 

6.  comments  were  liberally  added  to  make  the  program  more  user- 
friendly 
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V.  RESULTS  AND  DISCUSSION 


A.  SINGLE  TUBE 
I.  Repeatability 

Data  were  obtained  for  those  tubes  described  in  Chapter  III 
and  detailed  in  Tables  3.1  and  3.2.  Investigation  began  with  testing  of 
the  medium  diameter  tubes.  As  data  were  collected,  comparisons 
were  immediately  made  with  Van  Petten  [Ref.  6],  which  ensured 
validity  as  well  as  repeatability.  The  curves  on  all  figures  of  outside 
heat-transfer  coefficient  vs.  wall  temperature  difference  are  a  least 
-squares  fit  of: 

ho  =  aATr0,25  (5.1) 

Figures  5.1,  5.2,  and  5.3  compare  selected  data  for  medium  diameter 
tubes  under  vacuum  conditions  and  Figures  5.4,  5.5,  and  5.6  show 
similar  data  under  atmospheric  conditions. 

Figure  5.2  also  shows  typical  uncertainty'  bands  associated  with 
the  data.  Uncertainty  is  much  greater  for  the  low  cooling  water  flow¬ 
rates.  As  cooling  water  flow  rate  is  increased,  the  difference  between 
vapor  and  wall  temperature  increases,  the  outside  heat-transfer 
coefficient  decreases,  and  the  uncertainty  decreases.  Appendix  B 
gives  a  brief  decription  of  the  uncertainty  analysis  employed. 
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VAN  PETTEN 


VAN  PETTEN 


VAN  PETTEN 


VAN  PETTEN 


VAN  PETTEN 


VAN  PETTEN 


While  collecting  data,  any  sudden  variations  (from  expected 
trend)  in  the  outside  coefficient  could  have  been  indications  of  either 
dropwise  condensation  or  the  presence  of  noncondensables.  To  avoid 
dropwise  condensation,  cooling  water  flow  rate  was  initially 
maximized  so  that  condensate  fully  wetted  the  tube  surface  before  any 
data  were  taken.  This  method  seemed  to  be  critical  when  taking  data 
for  the  small  tubes,  and  especially  the  small  smooth  tube.  Several 
runs  were  conducted  for  the  small  smooth  tube  and  if  the  tube 
surface  had  not  been  initially  wetted  by  using  the  maximum  cooling 
water  rate,  'dry  patches  (potential  dropwise  sites)  on  the  tube 
surface  were  clearly  visible  and  the  heat-transfer  coefficients 
calculated  were  higher  than  expected.  Figures  5.7,  5.8,  and  5.9 
compare  selected  data  for  small  tubes  at  vacuum  and  Figures  5.10, 
5.11,  and  5.12  compare  small  tubes  at  atmospheric  pressure.  Van 
Petten  collected  data  for  the  QMCNPS  tube  only  at  atmospheric 
conditions.  Figure  5.13  is  a  comparison  of  that  data  with  that  of  this 
investigation. 

Table  5.1  summarizes  the  repeatabilty  of  data  with  that 
obtained  by  Van  Petten.  Comparisons  wen  based  on  the  difference 
between  the  least-squares-fit  calculaltions  of  the  heat  transfer 

coefficient  of  the  two  investigations  (%  difference  =  |(ho,  Coumes  -ho, 

Van  Petten)/ ho,Coumes]  x  100).  With  the  exception  of  tubes  2,  3  and  the 
QMCNPS  tube  at  atmospheric  conditions  (note,  all  are  small 
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Figure  5.10  Repeatability  of  Tube  3  at  Atmospheric  Condition 
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Figure  5.12  Repeatability  of  Tube  7  at  Atmospheric  Condition: 
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diameter  tubes),  the  data  are  in  good  agreement  The  discrepancies 
for  these  tubes  may  be  attributable  to  the  possible  presence  of 
dropwise  condensation  during  Van  Petten's  investigation.  As 
mentioned  previously,  small  diameter  tubes  were  more  susceptible  to 
the  dropwise  condensation  problem. 

2.  Inside  Heat -Transfer  Coefficient 

The  Sieder-Tate  type  coefficients  calculated  by  the  Modified- 
Wilson-Piot  technique  are  tabulated  in  Table  5.2.  Theoretically,  the 
value  of  Ci  for  tubes  with  the  same  internal  diameter  should  be  equal, 
unless  the  external  condensation  process  on  different  finned  tubes 
influences  the  circumferential  heat  flux  distribution.  However,  with 
the  exception  of  the  QMC  tube,  the  variation  observed  is  certainly 
within  the  combined  experimental  error  and  least-squares 
approximation  the  Modified-Wilson  Plot  technique.  The  relatively 
large  difference  exhibited  by  the  QMC  tube  is  probably  due  to  the 
insert.  Because  the  QMC  tube  has  a  larger  internal  diameter  (see 
Table  3.2),  the  insert,  which  was  designed  for  the  QMCNPS  tube,  did 
not  fit  as  tightly,  therfore,  the  enhancement  which  the  insert 
provided  was  not  as  great.  Table  5.3  shows  the  calculated  averages  of 
Ci  for  small  and  medium  diameter  tubes,  which  compare  favorably 
with  those  of  Van  Petten  [Ref.  6]. 
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Table  5.2 

SIEDER-TATE  COEFFICIENTS 


Tube 

Vacuum 

Atmospheric 

Small 

Smooth 

0.047 

0.048 

Finned 

2 

0.053 

0.050 

3 

0.053 

0.045 

4 

0.051 

0.047 

5 

(,).043 

0.044 

54* 

0.n46 

0.04? 

C 

0.045 

0.044 

7 

0.047 

0.04? 

QMC 

0.040 

0.036 

3MCNP5 

0.047 

0.046 

Medium 

Smooth 

0.07(! 

0.066 

Finned 

h 

0.066 

0.06s 

n 

0.067 

0.06  ' 

1 ;  , 

0.067 

0.06? 

;  • 

0.064 

0.065 

•  " 

0.065 

o.oc : 

K5 

0.064 

0.064 

*  Additional  da 

'.a  set  collected  (or 

Tube  5 

TABLE  5.? 

AVERAGE  VALUES  OF  C* 

Smalf  0.047  ±  0.003  (Van  Pettcn,  0.05!  i  0.00 

Medium  0.065  ±  0.002  (Van  Petten;  0.066  ±  0.00 


excludes  values  of  QMC  tube 


When  performing  calculations  to  determine  the  outside  heat- 
transier  coefficient,  the  average  values  of  Table  5.3  were  used.  The 
variation  is  small  and  if  the  Modified-Wilson  Plots  for  a  typical  tube 
are  compared,  there  is  very  little  difference  between  the  results 
obtained  by  using  the  values  of  Table  5.2  and  Table  5.3.  Figure  5.14 
shows  the  plots  of  a  randomly  selected  tube,  where  the  two  values  of 

Ci  have  been  used.  The  circles  are  representative  of  the  'best  fit 
value  (Table  5.2),  whereas  the  triangles  are  representative  of  the 
'average'  value  (Table  5.3).  The  X  and  lvalues  which  appear  as  axis 
values  are  those  which  appear  in  Equations  (4.13)  and  (4.12), 
respectively. 

When  working  with  instrumented  smooth  rubes,  Wanmaracchi 
et.  al,  [Ref.  71  found  that  data  were  better  correlated  when  a  constant 
term  (B)  was  added  tc  the  assumed  Sieder-Tate  type  relationship  for 
the  inside  heat-transfer  coefficient: 


Nu  =  CiRe0.2pr0,333(Mc/Mv)0.l4  +  B  (5.2) 

Equation  (5.2)  should  therefore  be  utilized  when  conducting  the 
Modified-Wilson  plot  technique.  It  would  entail  an  additional  iteration 
loop  (to  zero  in  on  B),  but  might  improve  the  data  fit  significantly. 
3.  Outside  Heat-Transfer  Coefficient 

The  outside  heat-transfer  coefficients  were  calculated,  as 
previously  stated,  using  Equation  (4.17).  Figures  5.15  and  5.16  depict 
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Figure  5.14  Modified  -  Wilson  -  Plot , 

'Best  Fit*  Versus  'Average'  Ci's  for  Tube  0  at  Vacuum 
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outside  heat-transfer  coefficient  data  versus  vapor-side  temperature 
drop  for  small  tubes  at  vacuum  and  at  atmosphere,  respectively. 
Figures  5.17  at  vacuum  conditions,  and  5.18  at  atmospheric 
conditions,  show  results  for  the  QMC  and  QMCNPS  tubes.  An  average 
Ci ,  from  Table  5.3,  was  used  to  obtain  the  QMCNPS  curves  while  the 
'best  fit'  value  cf  Table  5.2  was  used  for  the  QMC  curves.  The  large 
spread  in  the  data  can  probably  be  traced  to  the  determination  of  hi. 
Because  the  insert  did  not  appear  to  enhance  the  tube-side  of  the 
QMC  tube  as  efficiently  as  it  did  the  QMCNPS  tube  (the  larger  QMC  Di 
arguement),  the  QMC  hi  was  then  smaller,  therefore,  recalling 

Equation  (4.23),  its  h0  would  then  be  larger. 

Figures  5.19  and  5.20  are  the  vacuum  and  atmospheric  data  lor 
the  medium  diameter  tubes.  Also  plotted  on  all  the  figures,  and 
appearing  as  a  solid  line  below  the  tube  data,  is  the  Nusselt  theory 
prediction  for  a  smooth  tube. 
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Figure  5.16  Heat -Transfer  Coefficients  for  Small 
Tubes  at  Atmospheric  Conditions 
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Figure  5.17  Heat -Transfer  Coefficients  for  QMCNPS 
and  QMC  Tubes  at  Vacuum  Conditions 
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Figure  5.18  Heat-Transfer  Coefficients  for  QMCNPS 
and  QMC  Tubes  at  Atmospheric  Conditions 


85 


>  /K 


60 


O  Q  □  □  O  O 

ID  ^  CO  CM 

Ol 


Figure  5.19  Heat -Transfer  Coefficients  for  Medium 
Diameter  Tubes  Under  Vacuum  Conditions 


86 


wo 


60 


'  O  M  oo  CO 

f-H  «H  v-H  w  f-H 


0> 


I  II III 
i  mil 
mi  i 
mi  i 
mi  i 
mi  i 
i  mil 
i  mil 
I  II! !  ! 

+! 

+, 


I 

I 

I 

I 


d 

IP 

P 

/□ 


i  tin 

i  jury 

l*  * 
J*  4 

wi  i 
&r  i. 
ML  i 

d^4 

/ 

°  /  ////  / 

/  ////  / 

/  /^/  / 

/  jr/  / 

yyy  y 


i 

/□ 


o 

o 


o 

CO 

CO 

3 

2 


O 

in 


o 


o 

0) 


o 

(\J 


o 


Q 

(D 


Q 

in 


o 

't 

r\ 

O  . 

(T)  ll- 
I 

(I ) 
|— 

O 

CM 


-I  O 


O 


D 


Ol  /UM 


Figure  5.20  Heat  -  Transfer  Coefficients  for  Medium 
Diameter  Tubes  at  Atmospheric  Conditions 
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4.  Enhancement  Ratios 


Once  the  average  Ci  values  were  determined  all  data  were 

reprocessed  Using  an  average  Ci,  ex's  were  calculated  both  for  the 

smooth  tubes  (  <x*  )  and  for  the  finned  tubes  (  exf  ).  Enhancemer 
ratios  were  then  determined  using  Equations  (4.21)  and  (4.22). 

To  validate  the  accuracy  of  the  smooth  tube  data  it  was 
compared  to  Nusselt  theory,  and  also  to  the  empirical  correlation  of 
Fujii  et.  al  [Ref.  22].  Figures  5.21  and  5.22  show  small  smooth  tube 
data  at  vacuum  and  at  atmosphere,  respectively.  Figures  5.23  and  5.24 
show  medium  smooth  tube  data  for  vacuum  and  at  atmosphere, 
respectively.  The  smooth  tube  data  is  consistently  higher  than  that  of 
the  Fujii  correlation,  ranging  from  8  %  (medium  diameter  tubes)  to 
19  %  (small  diameter  tubes)  greater. 

Tables  5.4  and  5.5  summarize  the  enhancements  calculated 
as  well  as  the  retention  angle,  determined  by  Equation  (2.6),  and  the 
area  enhancement  The  area  enhancement  is  the  ratio  of  the  total 
surface  area  of  the  finned  tube  to  that  of  the  smooth  tube.  As  shown, 
heat-transfer  enhancement  is  greater  han  area  enhancement  with 
the  exception  of  the  closely  spaced  fins  (fully  flooded  conditions). 
However,  even  at  fully  flooded  conditions,  significant  heat-transfer 
enhancement  is  realized  Additionally,  atmospheric  conditions 
simultaneously  yield  greater  enhancements  and  smaller  retention 


s  *o.sy  /  nN 


Figure  5.24  The  Influence  of  Vapor  Velocity  on  Medium 
Smooth  Tube  Data  at  Atmospheric  Conditions 


92 


TABLE  5.4 


SUMMARY  OP  ENHANCEMENTS  AT  VACUUM  CONDITIONS 


Retention 

Angle 


Tubes 


Fin  Area 

Spacing  Enhancement 
(degrees)  (mm) 


Heat  Transfer 
Enhancement 


Small 


2 

180a 

0.25 

2.85 

1.85 

3 

180a 

0.50 

2.54 

1.71 

4 

156 

1.00 

2.15 

1.80 

5 

106 

1.50 

1.93 

1.64 

5A* 

106 

1.50 

1.93 

1.69 

6 

88 

2.00 

1.77 

1.75 

7 

59 

4.00 

1.46 

1.60 

Medium 

8 

180* 

0.25 

2.77 

1.83 

9 

180° 

0.50 

2.47 

1.79 

10 

1  10 

1.00 

2.10 

2.24 

i  i 

-A  X 

84 

1.50 

1.88 

2.39 

12 

71 

2.00 

1.74 

2.29 

13 

48 

4.00 

1.44 

2.16 

*  Additional  data  set  taken  for  Tube  5 

a  fully  flooded  condition  for  this  small  tube  occurs  at  a  fin 
spacing  of  0.96  mm 

b  fully  flooded  condition  for  this  medium  tube  occurs  at  a 
spacing  of  0.67  mm 
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TABLE  5.5 


SUMMARY  OF  ENHANCEMENTS  AT  ATMOSPHERIC  CONDITIONS 


Retention 

Angle 


Tubes 


Fin  Area  HeafTransfer 

Spacing  Enhancement  Enhancement 
(degrees)  (mm) 


Small 


2 

180c 

0.25 

2.85 

2.17 

3 

180c 

0.50 

2.54 

1.98 

4 

135 

1.00 

2.15 

2.27 

5 

98 

1.50 

1.93 

2.02 

5A* 

98 

1.50 

1.93 

2.01 

6 

82 

2.00 

1.77 

2.12 

7 

55 

4.00 

1.46 

1.94 

Medium 

8 

180d 

0.25 

2.77 

2.43 

9 

180d 

0.50 

2.47 

2.14 

10 

101 

1.00 

2.10 

2.50 

11 

78 

1.50 

1.88 

2.79 

12 

66 

2.00 

1.74 

2.47 

13 

45 

4.00 

1.44 

2.36 

'  Additional  data  set  taken  for  Tube  5 

c  fully  flooded  condition  for  this  small  tube  occurs  at  a  fin 
spacing  of  0.85  mm 

d  fully  flooded  condition  for  this  medium  tube  occurs  at  a 
spacing  of  0.60  mm 


angles  than  do  vacuum  conditions. 

Figures  5.25  and  5.26  are  plots  of  enhancement  ratio  versus 
fin  spacing  for  vacuum  and  atmospheric  conditions,  respectively. 
Small  diameter  tubes  become  fully  flooded  as  the  fin  spacing 
reduces  to  -0.9  mm,  whereas  medium  tubes  become  fully  flooded  at 
about  -0.6  mm.  Once  the  tube  is  fully  flooded,  enhancement  should 
increase  as  fin  spacing  is  decreased  further  because  metal  heat- 
transfer  surface  (with  a  igh  thermal  conductivity)  is  replacing  a  liquid 
film  (with  a  relatively  low  thermal  conductivity).  This  trend  may  be 
observed  for  the  medium  tubes  in  both  Figures  5.25  and  5.26. 
Consistent  with  Van  Petten  {Ref.  6],  Wanniarachchi  et  aL  [Ref.  7],  and 
Yau  et.  aL  [Ref.  8]  .  an  optimum  fin  spacing  of  1.50  mm  for  the  medium 
diameter  tubes  was  observed. 

The  small  tube  behavior  is  apparently  more  complicated. 
Considering  the  fully  flooded  region  (between  fin  spacings  of  0.25  - 
-0.9),  we  would  expect,  as  we  saw  for  the  medium  diameter  tubes,  to 
see  enhancement  increasing  as  fin  spacing  was  decreased.  However, 
what  we  observe  is  a  decrease  in  enhancement  between  fin  spacings  of 
1.0  -  0.5  mm  before  an  increase  is  seen  from  0.5  -  0.25  mm  The 
second  area  of  interest  is  the  data  at  a  fin  spacing  of  1.5  mm  As 
noted  in  Tables  5.4  and  5.5  an  additional  set  of  data  was  taken  for 
Tube  5.  This  was  an  effort  to  confirm  the  dip  which  the  small  tube 
enhancement  makes  at  a  fin  spacing  of  1.50  mm  This  dip  was  also 
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Figure  5.25  Heat  -  Transfer  Enhancement  Ratios  Under 

VacuunConditions 
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auLve.  wanaoHVBNa 

Figure  5.26  Heat -Transfer  Enhancement  Ratios  at 
AtmospherifConditions 
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recorded  by  Van  Petten.  When  an  uncertainty  of  about  10%  is  affixed 
to  the  small  tube  enhancement  ratios,  as  in  Figures  5.27  and  5.28,  an 
argument  could  be  made  that  the  dip  is  insignificant/nonexistent  and 
that  optimum  fin  spacing  for  the  small  tube  is  the  same  as  that  of 
medium  and  large  diameter  tubes,  namely  1.5  mm.  However,  Figures 
5.29  and  5.30  show  the  excellent  agreement  between  the  two 
data  sets  so  the  dip  may  be  real  Similar  dip  phenomena  have  been 
observed  for  R-113,  by  Marto  et  aL  [Ref.  42],  and  for  steam,  by  Yau  et. 
al.  [Ref.  8].  A  possible  solution  may  be  provided  by  the  three- 
dimensional  model  which  Honda  et.  aL  [Ref.  27]  have  proposed.  They 
theorize  that  the  dip  may  correspond  to  a  change  in  the  condensate 
film  profile  on  the  fin  root  tube  surface  in  the  unflooded  region  in 
going  from  large  spacings  to  small  spacings.  Further  investigation, 
both  experimental  and  theoretical,  is  therefore  required  to 
understand  this  small  tube  performance  characteristics. 

A  direct  comparison  of  the  enhancement  ratios  obtained  by 
Van  Petten  [Ref.  6]  can  be  somewhat  misleading.  Although  he 

performed  a  similar  reprocessing  of  data  he  used  average  Ci's,  to 
repiocess  and  calculate  <x's  for  his  finned  tubes,  but  did  not  use  those 

average  Ci's  to  calculate  new  tx's  for  his  smooth  tube.  Additionally, 

because  he  did  not  have  a  small  smooth  tube  he  used  the  medium 
smooth  tube  data  to  calculate  his  small  tube  enhancements.  The 
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Figure  5.27  Small  Tube  Enhancements  with  Uncertainties 
Under  Vacuum  Conditions 
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Figure  5.28  Small  Tube  Enhancements  with  Uncertainties 
at  Atmospheric  Conditions 
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Figure  5.30  Repeatability  ofTube  5  Data 
at  Atmospheric  Conditions 
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the  calculated  heat-transfer  enhancement  ratios  are  lower. 


TABLE  5.6 

SUMMARY  OF  SMOOTH  TUBE  tx'S 


Vacuum 


Atmospheric 


Small  Medium 
Van  Petten  0.803  0.803 

Coumes  0.933  0.838 

%  Difference  13.9  4.2 


Small  Medium 

0.867  0.867 

0.952  0.876 

8.9  1.0 


Theoretically,  as  Equation  (4.9)  shows,  the  effect  of  tube  diameter 


appears  only  in  the  value  of  F*>  and  oc  should  therefore  be  independent 
of  tube  diameter.  The  higher  value  obtained  for  the  small  tube,  both 
when  compared  to  the  medium  tube  data  and  to  the  Fujii  correlation 
could  be  attributed  to  undetected  dropwise  condensation  ,  or  perhaps 
simply  to  the  uncertainty  of  the  data. 

When  the  results  obtained  by  Van  Petten  [Ref.6]  for  large  tubes 
are  combined  with  results  presented  in  Figures  5.25  and  5.26,  it  can 
be  clearly  seen  that  the  enhancement  ratio  increases  as  tube  diameter 
increases,  as  shown  in  Figures  5.3 1  and  5.32. 
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5.  QMC  and  QMCNPS  Tubes 

Specific  discussion  of  these  results  has  been  previously 
deferred  One  of  the  thesis'  objectives  was  to  attempt  to  resolve 
differences  in  data  collected  during  a  concurrent  research  project 
being  conducted  by  Queen  Mary  College  (QMC)  of  the  University  of 
London. 

Marto  [Ref.  5]  compares  data  taken  at  NPS  for  R-113  and 
steam  with  the  data  taken  at  QMC.  The  NPS  data  was  gathered  for 
tubes  with  fin  heights  of  1  mm,  fin  thicknesses  of  1  mm,  and  a  fin  root 
diameter  of  12.7  mm.  while  the  QMC  data  were  for  tubes  with  a  fin 
height  of  1.59  mm,  fin  thicknesses  of  0.5  mm,  and  a  fin  root  radius  of 
12.6  mm.  As  expected,  and  shown  in  Figure  5.33,  because  of  the 
increased  heat-transfer  area  (due  to  larger  fin  height  and  fin  density), 
for  R-113,  the  QMC  data  showed  enhancements  greater  than  the  NPS 
data.  However,  for  steam,  just  the  opposite  trend  was  unexplicably 
observed. 

In  order  to  more  closely  correlate  the  data,  and  to  attempt  to 
resoh'e  the  conflicting  trends,  a  tube  was  manufactured  for  testing  at 
NPS  (termed  the  QMCNPS  tube).  Its  geometric  parameters  of  fin 
height,  thickness,  and  spacing,  internal  diameter  and  fin  root 
diameter  were  to  match  those  of  a  tube  being  used  in  the  QMC 
investigation.  Van  Petten  [Ref.  6]  reported  preliminary  results  for 
R-113  and  steam.  The  R-113  data  point  appears  on  Figure  5.33, 
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Figure  5.33  ComparisonofNPSandQMCDataforR-113 

|ReM2] 
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denoted  as  QMCNPS,  and  is  in  excellent  agreement  with  the  QMC  data 
of  Masuda  and  Rose  [Ref.  40).  However,  with  the  steam  data  point, 
there  was  no  agreement 

During  this  study,  a  tube  actually  used  in  the  QMC  study  was  also 
modified  and  then  tested  (termed  the  QMC  tube).  A  detailed 
geometric  comparison  between  that  tube  and  the  QMCNPS  tube  was 
previously  given  in  Table  3.2.  Figures  5.17  and  5.18  show 
comparisons  of  the  two  tubes'  data 

When  calculating  the  enhancement  ratios  for  the  tubes,  several 
methods  were  considered.  Table  5.2  lists  the  Sieder-Tate  coefficients 
for  all  tubes  tested.  Because  the  coefficients  for  the  QMC  tube  were 
considerably  lower  than  that  of  the  QMCNPS  tube  (and  the  other 
small  tubes),  two  different  enhancements  were  calculated  and  are 
reported  here.  One  ratio  was  based  on  the  'best  fit*  coefficient  of 
Table  5.2,  while  the  other  used  the  'average'  coefficient  of  Table  5.3. 
Figures  5.34  and  5.35  show  the  data  plotted  using  the  '  best  fit  '  and 
'average'  coefficients,  under  vacuum  and  at  atmospheric  conditions. 
As  can  be  observed  the  'average'  coefficient  with  the  least-squares-fit 
distorts  the  data  considerably,  when  compared  to  the  '  best  fit*  curve. 
Table  5.7  lists  the  enhancement  ratios  determined  using  both 


schemes. 
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Figure  5.34  Data  for  the  QMC  Tube  Under  Vacuum 
Conditions  Using  1  Best  Fit '  Versus  '  Average '  Values  of  Ci 
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Figure  5.35  Data  for  the  QMC  Tube  Under  Vacuum 
Conditions  Using  '  Best  Fit'  Versus  'Average'  Values  of  Ci 
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TABLE  5.7 


ENHANCEMENT  RATIOS  FOR  QMC  AND  QMCNPS  TUBES 

Vacuum  Atmospheric 

QMC  QMCNPS  QMC  QMCNPS 

Best  Fit  2.59  1.98  3.54  2.68 

Average  2.05  2.05  2.64  2.61 

Although  the  'average'  method  gives  enhancement  ratio  results  that 

are  in  close  agreement.  Figures  5.34  and  5.35  graphically 

demonstrate  the  reason  why  it  should  not  be  done.  Figure  5.36 

shows  the  enhancement  ratios  for  the  'b«.  fit*  scheme  at  atmospheric 

conditions,  also  included  are  actual  QMC  data  [Ref.  43]  for  steam. 

The  inconsistency  of  the  different  data  is  clearly  seen,  but  is  still  not 

fully  understood  Possible  causes  could  include: 

1.  discrepancies  in  geometry  of  tubes,  detailed  in  Table  3.2 
and  Figures  3.11  and  3.12 

2.  calculation  of  enhancements 

3.  use  cf  spiral  insert  and  its  effect  on  the  value  of  the  inside  heat- 
transier  coefficient 

The  last  item  appears  to  be  the  most  obvious  starting  pla  e  to  begin  a 
systematic  review  of  the  difference  in  the  two  investigations.  The 
QMC  data  for  steam  was  obtained  without  using  an  insert  Therfore, 

oc's  for  both  smooth  and  finned  tubes  could  be  significantly  different 
Also,  the  good  agreement  for  R-113  may  be  due  to  the  fact  that  the 
data  obtained  at  NPS  was  also  without  an  insert 
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Figure  5.36  '  Best  Fit '  Enhancement  Ratios  at  Atmospheric 

Conditions  for  QMC  Tubes  and  QMCNPS  Tube 
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B.  INUNDATION  EFFECT 


Flow  visualization  studies  of  the  effect  of  condensate  inundation 
were  conducted  for  both  smooth  and  finned  tubes.  Preliminary  static 
(i.e.,  non-condensing)  studies  with  water  had  demonstrated  that  the 
design  of  the  inundation  tubes,  Figures  3.6  and  3.7,  would  support  the 
inundation  rates  required  to  observe  the  full  range  of  flow  modes 
which  Honda  et.  aL  [Ref.  33]  had  described. 

The  condensate  return  system,  in  its  original  design  was  unable  to 
maintain  a  constant  inundation  rate,  because  of  pump  cavitation.  A 
modified  system  was  temporarily  installed  to  provide  recirculated 
condensate  to  a  single  inundation  tube  in  order  to  obtain  photographs 
for  the  initial  phase  of  flow-visualization.  Two  medium  smooth  tubes, 
one  in  the  dummy  tube  position,  the  other  in  the  active  tube  position, 
were  placed  in  the  test  section  and  atmospheric  operating  conditions 
were  reached  Visual  observations,  made  as  inundation  rate  was 
increased  revealed  the  four  falling  modes  .  A  video  camera  was  used 
to  record  the  condensate  flow  patterns  as  the  inundation  rate  was 
varied  Figures  5.37  through  5.42  are  Polaroid  snapshots  showing 
the  various  modes,  as  well  as  a  value  for  the  dimensionless  parameter, 
K,  described  in  Equation  (2.26).  In  the  droplet  and  column  modes, 
the  condensate  falls  at  regularly  spaced  intervals  (corresponding  to 
the  Taylor  wavelength,  Yung  et  aL  (Ref.  44]),  and  spreads  axially  over 
the  entire  surface  of  the  smooth  tube.  At  the  maximum  inundation 


113 


5.37  Inundation  of  Smooth  Tube, 
Droplet  Mode,  K  *  0.017 


Figure  5.38  Inundation  of  Smooth  Tube, 
Droplet  Mode,  K  =  0.045 
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Figure  5.39  Inundation  of  Smooth  Tube, 
Droplet  -  Column  Mode,  K  =  0.095 
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F igure  5.40  Inundation  of  Smooth  Tube, 
Column  Mode,  K  =  0.15 


Figure  5.4 1  Inundation  of  a  Smooth  Tube, 
Column  Mode,  K  =  0.21 
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Figure  5.42  Inundation  of  Smooth  Tube, 
Column  -  Sheet  Mode,  K  =  0.27 
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rate,  flow  nas  nearly  transitioned  completely  to  ihe  sheet  mode. 

The  smooth  tubes  were  removed  and  a  pair  of  medium  diameter 
Tinned  tubes  (  1  mm  fin  spacing,  thickness  and  height)  were  installed 
in  their  places.  Atmospheric  conditions  were  achieved  and  visual 
observations  were  again  made.  Both  35mm  and  video  cameras  were 
employed  to  record  the  failing  modes.  Figures  5.43  through  5.52  are 
35mm  prints  which  show  the  flow  pattern  and  the  corresponding 
.aiue  of  X.  Just  as  in  the  smooth  tube  case,  in  the  droplet  and  column 
modes  the  condensate  falls  at  regularly  spaced  intervals.  However,  as 
the  photographs  clearly  show  the  condensate  falls  to  the  tube  below 
unci  then  flews  circumferentially  down  the  grooves  between  several 
adjacent  fins.  This  is  in  contrast  to  the  smooth  tubes  where  the 
condensate  has  a  tendency  to  flow  axially  as  well  as  circumferentially 
around  rhe  tube.  The  space  between  these  condensate  inundated 
regions  is  unaffected,  and  therefore  a  significantly  enhanced  heat- 
.ranoier  performance  is  expected  when  compared  to  the  smooth  tube, 
whose  entire  surface  is  wetted  by  the  axial  migration  of  rhe  falling 
condensate. 

Honda  et.  aL  [Ref.  33]  reported  K  values  for  failing  mode 
transitions.  A  comparison  of  those  reported  values  (for  low  surface 
tension  fluids)  with  those  observed  during  this  preliminary 
investigation  are  recorded  in  Table  5.8. 
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Figure  5.43  Inundation  of  Finned  Tube 
Droplet  Mode.  K  =  0.030 
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f  igure  5.44  Inundation  of  Finned  Tube, 
Droplet  Mode,  K  =  0.065 
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Figure  5.45  Inundation  of  Finned  Tube, 
Droplet  -  Column  Mode,  K  =  0. 100 
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Figure  5.46  Inundation  of  Fnned  Tube, 
Column  Mode,  K  =  0. 14 


Figure  5.47  Inundation  of  Finned  Tube 
Column  Mode,  K=  0.17 
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Figure  5.48  Inundation  of  Finned  Tube 
Column  Mode.  K  =  0.20 


Figure  5.49  Inundation  of  Finned  Tube 
Column  Mode,  K  =  0.24 
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Figure  5.50  Inundation  of  Finned  Tube 
Column  Mode,  K  =  0.27 


128 


Figure  5.5 1  Inundation  of  Finned  Tube 
Column  Mode,  K  =  0.31 
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Figure  5.52  Inundation  of  Finned  Tube 
Column  -  Sheet  Mode,  K  =  0.34 
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TABLE  5.8 


DIMENSIONLESS  PARAMETER.  K 


Honda  Smooth 
Tubes 


Falling  Modes 


Droplet  -  Droplet/Column 
Droplet/Column  -  Column 
Column  -  Column/  sheet 


0.06-0.13  0.045 

0.06-0.13  0.095 

0.37-0.47  0.27 


Finned 

Tubes 


0.065 

0.10 

0.34 


VL  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

1.  The  heat  transfer  enhancement  of  integral-fin  tubes  is  affected 
by  root  diameter.  Enhancement  increases  as  root  diameter  increases. 

2.  The  optimum  fin  spacing  is  near  1.5  mm  for  the  medium  tube, 
as  found  in  previous  published  investigations. 

3.  Small  tube  behavior  deviates  from  the  standard  optimum  fin 
spacing,  and  appears  to  exhibit  an  optimum  fin  spacing  of  1  mm. 

4.  Data  taken  for  the  QMC  tube  is  unexplainably  higher  than  that 
obtained  by  the  joint  project  at  Queen  Mary  College.  The  utilization  of 
a  spiral  insert  during  the  tests  at  NPS  may  be  the  cause  . 

5.  Condensate  inundation  flow  modes  were  correlated  to  the 
dimensionless  parameter  K. 

6.  Photographic  and  visual  observations  clearly  demonstrated  axial 
migration  of  condensate  for  smooth  tubes,  resulting  in  a  totally  wetted 
surface.  For  finned  tubes,  circumferential  migration  (limited  to  several 
channels  between  adjacent  fins)  was  observed,  leaving  unaffected 
portions  of  the  finned  tube  free  of  inundation. 

B.  RECOMMENDATIONS 

1.  Obtain  small  tubes  with  fin  spacings  of  1.25  and  1.75  mm  in 
order  to  determine  an  optimum  fin  spacing  for  this  small  diameter 
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tube. 


2.  Re-examine  all  the  differences  between  the  tubes  used  in  the 
joint  program  with  Queen  Mary  College  and  then  repeat 
measurements  in  an  effort  to  clear  up  the  inconsistencies  in 
the  steam  data 

3.  Take  data  with  and  without  the  spiral  insert  and  process  the 
resulting  data  with  the  new  form  of  the  Sieder-Tate  relationship 

(i.e.  B»*0). 

4.  Replace  the  mercury-in-glass  manometer  with  direct  reading 
pressure  gage. 

5.  Relocate  the  pressure  and  temperature  probes  in  the  inundation 
test  section  so  they  are  upstream  of  the  test  tubes  instead  of 
downstream  as  presently  located 

6.  Obtain  additional  finned  tubes  to  support  further  condensate 
inundation  studies. 

7.  Conduct  condensation  studies  with  condensate  inundation 
effects  for  various  fluids,  including  steam  and  R-l  13. 
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APPENDIX  A 


SYSTEM  START-UP  AND  SHUT-DOWN  PROCEDURES 

The  system  is  started  in  the  following  manner 

1.  Fill  the  boiler  with  distilled  water  to  a  level  of  four  to  six  inches 
above  the  top  of  the  heating  elements. 

2.  Check  purge  pump  system;  oil  level,  drain  valves,  system  vent 
valve. 

3.  Turn  on  data  acquistion  program  and  load  software  program 
entitled  DRP12S. 

4.  Open  fill  valve  for  cooling  water  sump  supply. 

5.  Turn  on  ooling  water  supply  pumps  and  adjust  tube  flow  rate  to 
approximately  45%. 

6.  Open  valve  from  tap  water  system  to  auxiliary  condenser  and 
adjust  flow  rate  to  approximately  40%. 

7.  Check  system  for  water  leaks. 

8.  Energize  heaters  and  adjust  voltage  to  -50  volts  for  system 
warm-up. 

9.  If  conducting  vacuum  run,  energize  vacuum  purge  pump,  then 
open  system  vent  valve  (single  turn). 

10.  Monitor  sytem  temperature  (via  thermocouple  voltage  reading) 
and  pressure  (mercury-in-glass  manometer)  carefully. 

11.  If  conducting  vacuum  run  wait  until  system  pressure  reaches 
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~120  mm  Hg  (left  hand  column)  and  then  gradually  increase  voltage  in 
10  volt  increments  to  90  volts.  Obtain  desired  operating  conditon  by 
manually  controlling  flow  through  auxiliary  condenser. 

12.  If  conducting  atmospheric  run,  open  system  vent  valve  and  purge 
system  drain  valve.  After  - 1/2  hour  warm-up  begin  increasing  voltage 
in  10  volt  increments  until  ~95  volts  is  reached.  Wait  until  system 
pressure  drops  (just  as  hot  vapor  bubble  makes  its  180°  turn  at  the  top 
of  the  apparatus),  then  close  system  vent  valve  and  purge  system 
drain  valve.  Energize  vacuum  pump  and  open  the  system  vent  valve 
~  1/3  turn.  Simultaneously  raise  voltage  in  20  volt  increments  to  175 
volts  while  throttling  flow  through  the  auxiliary  condenser  to  obtain 
desired  operating  conditions. 

13.  Adjust  cooling  water  flow  meter  setting  to  proper  %  for  initial 
data  run. 

14.  Allow  the  vacuum  pump  to  operate  continuously,  monitor  oil 
level  through  sight-glass. 

The  system  is  secured  in  the  following  manner 

1.  Secure  power  to  the  boiler  elements. 

2.  Isolate  and  secure  the  purge  system. 

3.  Slowly  open  the  system  vent  valve  and  purge  drain  valve  in  order 
to  allow  the  system  to  come  to  atmospheric  pressure. 
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4.  Continue  to  circulate  cooling  water  through  the  tube  and 
auxiliary  condenser  to  assist  in  cooling  down  the  system. 

5.  Monitor  system  temperature,  and  when  cool  turn  off  data 
aquisition  system 

5.  Secure  the  cooling  water  systems. 

6.  Check  oil  level  (and  quality)  in  purge  pump. 

7.  Periodically  change  distilled  water  in  boiler. 


136 


APPENDIX  B 


UNCERTAINTIES  AND  ASSUMPTIONS 

Uncertainties  are  unavoidably  introduced  whenever  taking 
measurements.  They  are  dependent  upon  the  accuracy  of  measuring 
equipment  and  the  operator's  skill  in  using  them.  Although 
uncertainties  associated  with  single  measurements  may  be  small, 
when  combined  with  others  in  the  data  reduction  scheme,  they  can 
quickly  propagate  into  a  rather  large  uncertainty.  The  uncertainties  of 
the  measured  quantities  in  this  investigation  were  determined  by 
using  the  follow  relationship  proposed  by  Kline  and  McClintock  [Ref. 

ir\ 

-to,. 

Wr  =  {  |(6R./oXl)WlJ-  +  |(6R/6X2)W2]£  +  -  +  [(6R/6Xn)Wrt|£ 

(B.l) 

where 

R  =  the  result  of  the  calculation, 

Wi  =  the  uncertainty  of  the  result, 

Xa  =  the  measured  independent  variables  (i=  1,2,3 . n). 

Wi  =  the  uncertainty  in  the  measured  variable  (i=  1,2,3 . n). 

Georgiadis  [Ref.  35]  provides  a  detailed  account  of  the  development  of 
the  uncertainty  analysis.  Mitrou  [Ref.  46)  created  a  program  which 

calculated  uncertainties.  Sample  outputs  follow; 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name: 

25V 

Pressure  Condition: 

Vacuum 

(  1  1 

tPa  ) 

Vapor  Temperature 

= 

48.40 

( Deg  C  ) 

Water  Flow  Rate  ( %  ) 

= 

GS.00 

Water  Velocity 

= 

G  .49 

( m/s  ) 

Heat  Flux 

= 

2 . 359E+05 

( W/mA2  ) 

Tube-metal  thermal  conduc. 

= 

385.0 

(W/m.K ) 

Sieder — Tate  constant 

= 

0.0470 

UNCERTAINTY  ANALYSIS: 

VARIABLE 

PERCENT  UNCERTAINTY 

Mass  Flow  Rate  ,  Md 

0.95 

Reynolds  Number,  Re 

1  .28 

Heat  Flux  ,  q 

1.28 

Log-Mean-Tem  Diff,  LMTD 

.73 

Wall  Resistance,  Rw 

2.G3 

Overall  H.T.C.  ,  Uo 

1  .47 

Water — Side  H.T.C.  ,  Hi 

G .  48 

Vapor — Side  H.T.C.,  Ho 

B .  9 1 

OATA  FOR  THE  UNCERTAINTY  ANALYSIS: 

File  Name:  12SA 

Pressure  Condition:  Atmospheric  (101  kPa  ) 

Vapor  Temperature 

=  100.05 

( Deg  C  ) 

Water  Flow  Rate  <%) 

=  54.00 

Water  Velocity 

=  2.99 

( m/s  > 

Heat  Flux 

=  8 . 792E+05 

<  W/m'  2  ) 

Tube-metal  thermal  conduc. 

=  385.0 

(W/m.K  ) 

Sieder-Tate  constant 

-  0.0G50 

UNCERTAINTY  ANALYSIS: 

VARIABLE  ’ERCENT  UNCERTAINTY 


Mass  Flow  Rate ,  Md 

1  .16 

Reynolds  Number,  Re 

1  .42 

Heat  Flux,  q 

1 .25 

Log-Mean-Tem  Diff,  LMTD 

.16 

Wall  Resistance,  Rw 

2.G7 

Overal 1  H.T.C.  ,  Uo 

1  .27 

Water-Side  H.T.C.  ,  Hi 

3.30 

Vapor-Side  H.T.C.,  Ho 

5.98 
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Perhaps  just  as  important  as  uncertainties  in  the  final  evaluation 
of  collected  data  is  considering  the  assumptions  used  throughout  the 
data  reduction  analysis.  A  partial  list  of  the  assumptions  used  during 
this  investigation  follows: 

•  no  partial  dropwise  condensation  occurs 

•  fouling  not  included  as  thermal  resistance 

•  noncondensable  gases  not  included  as  thermal  resistance 

•  strictly  radial  conduction  through  the  tube  wall(i.e.  concentric 
isotherms  were  assumed) 

•  no  circumferential  variation  in  wall  temperature  profile 
considered 

•  thermal  properties  of  tube  material  considered  uniform  and 
constant 

•  thermal  resistance  of  black  oxide  layer  considered  negligible 

•  modelling  of  inlet  and  outlet  portions  of  tube  as  fins 

•  evaluation  of  cooling  water  properties  at  bulk  average 
temperature 

•  use  of  Modified-Wilson-Plot  technique 

•  assumption  of  Sieder-Tate  type  relationship  for  inside 
heat-transfer  coefficient  with  swirl-type  insert 

•  assumption  of  Nusselt-type  relationship  for  outside  heat- 
transfer  coefficient 
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•  calcualtion  of  condensate  film  temperature  as 
(2/3)Tw*ll  +  (l/3)T**t 


140 


APPENDIX  C 


RAW  DATA 

Table  C.l  contains  the  names  of  the  raw  data  files  with  the 


corresponding  tube  designation.  A  sample  file  follows  the  table. 

TABLE  C.l 

SUMMARY  OF  RAW  DATA  FILES 
SINGLE  TUBE  STUDY 

Tube  Vacuum  Atmospheric 


SmallTubcs 

Smooth 

SMV 

SMA 

2 

2SV 

2SA 

3 

3SV 

3SA 

4 

4SV 

4SA 

5 

5SV 

5SA 

5A 

5ASV 

5ASA 

6 

6SV 

6SA 

7 

7SV 

7SA 

QMCNPS 

QNPSV 

QNA 

QMC 

QMCV 

QMCA 

MediumTubes 

Smooth 

MEDV 

MEDA 

8 

8SV 

8SA 

9 

9SV 

9SA 

10 

10SV 

10SA 

11 

11SV 

11SA 

12 

12SV 

12SA 

13 

13SV 

13SA 
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Tube  Number:  08 

File  Name:  8SA 

Pressure  Condition:  Atmospheric 
Steam  Velocity:  1.0  (m/s) 


Data 

Vw 

Tin 

Tout 

Ts 

# 

<  m/s  ) 

(C) 

<C  ) 

(C  ) 

1 

1.16 

23.  17 

30.77 

99. 

91 

2 

1.16 

23.22 

30.87 

99. 

94 

3 

1  .48 

23.17 

29.67 

99. 

98 

4 

1.48 

23.21 

29.70 

99. 

93 

S 

1  .97 

23.09 

28.62 

99. 

97 

6 

1.97 

23.05 

28.58 

99. 

90 

7 

2.51 

22.97 

27.76 

99. 

96 

8 

2.51 

22.97 

27.78 

100. 

22 

9 

2.99 

22.89 

27.23 

99. 

92 

10 

2.99 

22.90 

27.24 

99. 

92 

1  1 

3.42 

22.84 

26.89 

99. 

88 

12 

3.42 

22.85 

26.99 

99. 

94 

13 

3.85 

22.80 

26.56 

99. 

92 

14 

3.85 

22.77 

26.55 

100. 

06 

15 

4.39 

22.56 

26.00 

100. 

00 

16 

4.39 

22.52 

25.96 

99. 

83 

17 

1.16 

23.20 

30.91 

100. 

03 

18 

1.16 

23.22 

30.94 

100. 

14 
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